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PREFACE 




To say that there have been many articles written over the past 
few years dealing with phase-locked devices seems almost an under- 
statement. Besides this work, I am aware of two books being written 
on the subject, those of Viterbi and Van Trees, and there may be still 
others. This monograph was prompted by a need at the Jet Propulsion 
Laboratory for a practical handbook for the design, testing, evalua- 
tion, and philosophy of phase-locked systems. Such a book should 
possess several qualities: it should contain all the formulas needed for 
the accurate design of a complicated receiver; it should provide a 
unified notation with unambiguous definitions of receiver parameters; 
it should tell how to pick these parameters, based on any specific job 
required of a communications system; it should give special attention 
to the general system design problems and procedures, with all the 
details included; it should serve as a reference against which a receiver 
in the field can be checked to see whether or not it is working prop- 
erly, and thereby meet its proper specification; and, finally, it should 
contain that philosophy which has evolved to ensure the successful 
fabrication of the most sensitive, flexible, and stable receiver in the 
world today. It has been my intention to fulfill these goals in the pages 
you see here. 

The approach I have followed in trying to assemble words and 
formulas enough to succeed in my aim is that the presentation should 
expose the reader to the reasoning by which the theory evolved and 
to a few of the steps necessary for him to follow this evolution mathe- 
matically, and should set forth its results in as concise a treatment as 
precision could permit. 

Thus, not only in intent, but also in approach, is the treatment in 
the following pages different from that to be found elsewhere. The 
two books mentioned above are theoretical treatments, rigorously 
argued, and excellent reading for one desiring specific insight to the 
purer aspects of communications. But the content of these works 
necessarily ends before most of the detailed systems engineering tech- 
niques are evoked. For example, the band-pass-limiter loop is no- 
where treated with the detail it is here. 


XI 





»w <vatt?n t - ^ T“ ? ^ tn-^au 


JPL TECHNICAL REPORT NO. 32-819 

PREFACE (Cont’d ) 

While this is intended as a working book, by no means have I 
presumed to consider the practical aspects of building a better VCO, 

RF mixer, IF amplifier, etc. These are components that develop nat- 
urally with the state of the communications art. What is given is / 

meant to provide the systems engineer with the tools he needs to effect 
a successful analysis of his particular job. 

Any exposition of a broad topic mathematically is apt to frustrate 
the conscientious author almost beyond comprehension in the matter 
of notation. He soon exhausts both Roman and Greek alphabets, and, 
as the frustration becomes extreme, he begins to eye script, Russian, 

Hebrew, and Germanic characters in the desperate hope that he can 
use separate symbols for all the quantities to be defined, before all of 
these, too, are used up. Suffice it to say that I have retained only the 
Roman and Greek at some awkwardness, being forced away, in some 
cases, from a common notation because of multiple demands on the 
same letter. An attempt was made to resolve these conflicts by giving 
the more important quantity precedence. I have also tried to keep 
subscripts as simple as possible by keeping subscripted subscripts and 
other such complexities to a minimum. 

The subject is treated in two volumes, the first of which consists of 
ten chapters dealing primarily with the detailed tracking aspects 
of phase-locked devices. Topics involving modulation in various modes 
are relegated to the second volume. Each of the volumes is further 
divided into two portions: the first, a section setting forth theory and de- 
sign optimization; and the latter, a summary of formulas and methods. 

It is hoped that the reader will find in these pages as much insight into 
the mysterious behavior of the phase-locked loop as the author sup- 
poses he developed in writing them. 
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ABSTRACT 

This is Volume I of a two-volume work on the theory of phase- 
locked receivers, with pertinent reference material on practical re- 
ceiver design. Volume i is primarily devoted to the performance of 
carrier-tracking loops, including a rigorous treatment of narrow-band 
systems having IF limiters. The bulk of the work is based on a theo- 
retical linear model, but a nonlinear method is also presented to predict 
behavior near the threshold. The emphasis throughout the work is 
toward completeness, simplicity, and internal consistency of the ma- 
terial assembled. 

Part I of this Volume is an exposition of the theory, the resulting 
equations, and the design philosophy that have enabled the phase- 
lock concept to evolve into the basic principle underlying the most 
sensitive receivers in the world today. Part II is a condensed version of 
Part I, intended as a quick reference to formulas, definitions, and 
salient design considerations. 
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PART I 

THE THEORY OF PHASE-LOCKED LOOPS 
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CHAPTER 1 

INTRODUCTION AND HISTORY OF THE PHASE-LOCKED LOOP 


While the origins of automatic phase control date back 
to the 1920’s and 30’s, the first serious application of the 
concept began as a horizontal-line synchronizing device 
for television in the late 40’s. Shortly thereafter, Jaffee and 
Rechtin showed how a phase-locked loop could be used as 
a tracking filter for a missile beacon, and how the loop 
parameters could best be specified. The first analysis 
including the effects of noise appeared in a paper pub- 
lished by Jaffee and Rechtin in 1955. 

Basically, a phase-locked loop is an electronic servo- 
mechanism that operates as a coherent detector by con- 
tinuously correcting the frequency of its local oscillator 
according to a measurement of the error between the 
phase of the incoming signal and that of its local oscillator. 
The simplest form of loop is shown in Fig. 1-1. The precise 
relationship between the input and response functions is 
a nonlinear integro-differential equation from which very 
little information concerning loop behavior is analytically 
available, in the general case at least. With a very few 
nonrestrictive assumptions mathematically the configu- 
ration given in Fig. 1-1 can be replaced by that in Fig. 1-2 
(we shall indicate precisely why this is so a little later). 
This model first appeared, without proof, in a paper by 
Develet in 1956, but in the absence of noise it had been 
used by several authors, notably Gruen. 

Viterbi found solutions for a number of loop filters and 
various input frequency functions in the no-noise case by 
analog simulation. 



Fig. 1-1. Basic configuration of a simplo phast-lockod 
loop. The mixor output, filtered by F (•), is usod 
to control tho frequency of the voltage- 
controlled oscillator (VCO). 


n[f ) - 2 l/z /i,(/)cos {uit +§) 



Fig. 1-2. A mathematically equivalent model of the 
simple phase-locked loop. Here there are two 
separate inputs, 6 and nit), whereas the two 
are combined in the actual input. 

The case in which additive noise is present has been 
treated by a variety of approximate methods. The first 
approach, by Jaffe and Rechtin, essentially replaced the 
sinusoidal nonlinearity of the model of Fig. 1-2 by a linear- 
amplifier of gain A, a case applicable when the phase 
error is very small. Margolis analyzed the nonlinear oper- 
ation in the presence of noise by perturbation methods, 
obtaining a series solution for the loop differential equa- 
tion, and, using only the first few terms of the series, he 
determined approximate moments of the phase error. 
Develet applied Booton’s quasi-linearization technique to 
replace the sinusoidal nonlinearity by a linear amplifier 
whose gain is the expected gain of the device. More 
recently. Van Trees obtained a Volterra series representa- 
tion of the closed-loop response by a perturbation method 
similar to the method employed by Margolis, but with the 
advantage of the simplified model he obtained more exten- 
sive results. Fokker-Planck, or continuous random-walk, 
techniques yield exact expressions for the statistics of the 
random phase error process. Unfortunately, expressions 
in closed form are available only for the first-orde? loop 
(i.e., when the filter is omitted). Such techniques were first 
applied to this problem by Tikhonov, who was able to 
determine the steady-state probability distribution of the 
first-order loop phase-error and an approximate expression 
for the distribution when the loop contains a one-stage 
RC filter. Viterbi extended this work on the first-order 
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loop also to obtain the mean time to loss-or-lock (giving 
the frequency of skipping cycles). 

In what follows, we shall review many of these analyses 
and add another, due to the author, wherein the spectral 
density of the phase process is approximated. This method, 
incidentally, is conceptually as simple as the linear approx- 
imation methods and yields results almost as exact as the 
results gotten by Fokker-Planck techniques over the entire 
useful range of the device. 

These analyses have not been merely of academic 
interest. They have paved the way for building the most 


narrow-band, sensitive, flexible receivers in the world. 
Phase-locked loops are used as filters to “clean up” the 
output of frequency multipliers. The phase-lock principle 
has been used in ranging devices and in radar systems 
capable of tracking a planet (Venus) with range-jitter less 
than 500 meters! It has been used to synchronize telemetry 
data, to derive bit and word synchronization . . . and the 
list could go on at some length. 

Because there is such a long list of jobs it can do so well, 
it is only natural that the phase-locked loop has received 
attention. Significant portions of the fruits of this attention 
are revealed in the ensuing chapters. 
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CHAPTER 2 

FUNDAMENTAL CONCEPTS 


In this Chapter we shall define such concepts as 
auto- and cross-correlation functions, spectral density, 
bandwidths, etc. This defining must be done carefully, 
as it can lead to one of the prime sources of reader con- 
fusion. We shall assume that the reader is familiar with 
certain statistical properties of signals and noise, so that 
the exposition here need not be a tutorial one, but merely 
descriptive. 


2-4. Statistics 

First, consider a single time function x(t). Ask yourself: 
Can x(t) be specified exactly at every value of t? If it is, 
x(t) is said to be deterministic. Otherwise, there is un- 
certainty about the character of x(t). In such cases, any 
partic alar observed waveform x(t) is said to be a sample 
function of a random process (others call x(t) a stochastic 
process). There are varying degrees of randomness one 
can think of here. For example, x(t) = ccs t is a well 
defined, deterministic process. However, x(t) — A cost, in 
which A can assume any arbitrary value (according to 
some probability law) is a sample function from a random 
process, even though we can measure x(0) = A, and from 
then on, x(t) is known exactly. The randomness here is 
evident, for if we were presented with another sample 
function from the same process, we would probably have a 
different A. 

Another example that occurs frequently is the sample 
function x(t) = cos (t + $), where 6 can assume any value 
over [0,2>r] with equal likelihood. Again, 6 can be found for 
any particular observed x(t) very easily, and forever there- 
after x(t) is known exactly. However, there is a great differ- 
ence in the types of processes represented by A cos t and 
cos(t + 6). If one were to ask “what is the mean value of 
x(t) at time t = #,?” in the first case one would answer 
"E(A) cos #i," where E(A) represents the expected value 
of A with respect to its random elements: 

E(A) = J A p(A)dA (2-1) 

whereas, in the second, one would answer “zero," because 
for every value of t u 

E[*(*i)] = ^ J ” cos (t t +t)d$ = 0. (2-2) 


The averaging operator E here, we note, is not the same 
as the operator that averages over time. In fact, in the ( 
first case above, the mean value depends on f, and in the 
second, it does not. 

Generally, a process whose statistical behavior is inde- 
pendent of the time origin — as is cos(t+0) — is said to be 
a stationary process, while others not exhibiting this prop- 
erty are nonstationary. 

Another point that can be made is that, for A cos t, 
averages over the entire set, or ensemble, of sample func- 
tions are not the same as averaging with respect to time. 

A process in which the time averages involving all 
functions of x(t) are the same as averages taken over the 
ensemble of sample functions is said to be ergodic. This is 
a behavior different from stationarity, as witnessed by the 
fact that the function A cos(t+6) has its time-mean-square 
value equal to A 2 / 2, whereas its ensemble mean-square 
value is E(A 2 )/ 2. Ergodic processes are always 1 station- 
ary, but, as our example shows, the converse need not 
be true. 

If x(t ) is a well-behaved sample function, the expression 

- 1 im±J\(t)x(t+r)dt (2-3) 

defines the time-autocorrelation function of x(t). The 
entire process has a statistical analog 

= E [x(t,)x(t,)] (2-4) 

in which averaging is performed with respect to the 
statistical variables in x(t). When x(t) is from a stationary 
process, the latter is a function of t a — 1 2 , rather than ti and 
f, separately, in which case we write 

H f ,(r) = E[x(t)x(t+r)]. (2-5) 

Any process x(i), stationary or not, satisfying this last 
particular equation is called widesense stationary. Under 
the further restriction that x(t) be ergodic, we have 

R~(r) = £R„(t). (2-6) 


*Thii is not strictly true under a more formalized math em ati c al 
definition of ergodidty, but this formalism is of no concern here. 
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The spectral density o£ a function x(t) is defined as the 
Fourier transform of its time autocorrelation function, 


is the cross correlation function of the two processes; when 
stationarity prevails, we set — t: 


/ +oo 

(2-7) 

and the spectral density of a wide-sense stationary process 
is defined as 


/ +oo 

R ri (r)e-)^dt. (2-8) 

When x(t) is a function from a stationary ergodic process, 
these spectral densities are equal with probability one, 

S„(1<o) = S„(jw). (2-9) 

We shall usually not make such distinctions in notation 
in the future. For the most part, we deal with stationary 
ergodic processes whenever analysis requires a statistical 
treatment. 

Having found S„(/ u ), we may replace /<u by the complex 
frequency variable s; the resulting function S zz (s) is the 
analytic continuation of S„(f< «) to the entire complex plane. 
When S zz is known, R zz can be found by inverse trans- 
formation. It is usual to insert <» = 2nf in S„ and treat 
spectra as functions of f. (Sometimes care must be taken 
in this step; e.g., |u>| -*■ (—s 2 )'*.) 

The mean of x(t) (i.e., its expected behavior) is denoted 

p(t) = E[x(t)] (2-10) 

and its variance (i.e., the mean-square variation about 

/ i(t )), is denoted 

o>(t) = E{[x(t)- fi (t)y) 

= £[**(«)] -At) (Ml) 

For stationary processes, both a- and p are independent 
of t, and, if R*,(oo) exists, 

pI = R*»(oo) 

<r* ~ fl**(0) ~ K**(oo) (2-12) 

Equations (2-3) through (2-12) can be altered to yield 
functions of average products of two different processes, 
say x(t) and y(t). Then, for example, 

JU(M>) = £ [*(*i)y(*»)] 


R„(r ) = E[x(t)y(t-r)] 
= R vz (~t). 


i'he Fourier transform S^,(/w) of R zy {r) is then called a 
cross-spectral density. 


2-B. Linear Filtering 

Consider a device, which we shall label H as depicted 
in Fig. 2-1, having an input x(t), giving rise to an output 
y(t). Suppose first that the input is a very sharp pulse 
having unit weight; that is, an impulse function fa(t): 

(i forO^t^T 

8r(*)=j r (2-13) 

( 0 for all other values of t. 


As T approaches zero, 8 T (t) becomes infinite at the origin, 
and the result is a Dirac delta function 8(f). One must be 
very careful in treating functions of this sort mathemati- 
cally, because of the lack of uniform convergence. How- 
ever, for our applications it has been shown that in all 
practical calculations we may treat 8(f) as the limit of 
8r(t) as T— >0. 


*(f> 


H 


yif) 


Fig. 2-1. Filtering device 


The response of the device H to 8(f) is called the unit- 
impulse response of the filter, denoted h(t). That is, h(t) is 
the output (Fig. 2-2) of H at time t when an impulse was 
applied at time t = 0. The filter is called realizable if h(t) 
= 0 when t < 0 (i.e., no response before an input). The 
filter is said to be linear if superposition holds: 

y(t) = J x(t x ) h(t—t x ) dti 

= [ ~x(t-t l )h(t l )dt l . ( 2 - 14 ) 


WJ 
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* (/) 



t 



Fig. 2-2. Response of a realizable filter to a 
unit-impulse function 


Thi r law states that the output y at time t is the super- 
pcjition of the input t, seconds ago, weighted by the 
decay response of the filter after t, seconds. 

When x(t ) and h(t) have Fourier or Laplace transforms 
(call them X(s) and //($)), then for linear filters, y(t) has a 
transform Y(s) related to the others by 

Y(s) = X(s)H(s). (2-15) 

However, many random processes (notably noise) do 
not possess transforms, so one cannot perform a system 
analysis based on X(s). Instead, one can compute statistics 
such as the correlation function of the output, 

IWr) = j* l*Ht l )h{t i )R„(T+U-U)dUdt i 

(2-16) 

and this expression can be transformed to yield the output 
power spectral density of the process: 

S„(M = 

= |H(M|*S^M. (2-17) 

This equation reveals one of the most important facts 
about linear filters: the spectral density of the output 
process is merely the product of the input density times 
the filter's response. By analytic continuation, 

SmM = H(s)H(-s)S^s). (2-18) 


The power 2 appearing at the filter output can be found 
by inverse transformation 

P v = RJO) = ± f~ |H(/iu)| 2 SMd». (2-19) 

2-C. Noise Bandwidth 

Suppose n(t) is a zero-mean, ergodic random process of 
such a nature that n(t) and n(f+r) are completely un- 
correlated for every t 0; that is, = 0 for every 

r 0. Such a process is often called a white noise process. 
It is often the character of such noises that 

R„„(r) = AU(r). (2-20) 

The spectral density of such a noise is thus uniform in 
the frequency domain 

S nn (ju>) = N„ for all — oo < w < oo . (2-21) 

It should be pointed out that R„„(0) indicates n(t) has 
infinite power, while S »«(/«■>) shows that over any finite 
frequency range, the power is finite. The N 0 above can be 
thought of either as the weight of the correlation function 
impulse, or else as the uniform height of the spectral 
density of the noise. 

Now consider what happens when a white noise n(t) is 
put into a linear filter H(s): the output spectral density is 

S„(M = N„|H(M| 2 (2-22) 

and the output noise power N from the filter is 

N = JjH(M| 2 d»]. (2-23) 

Whenever the integral in brackets converges, the output 
noise power is brought down from an infinite value to a 
finite one by the action of the filter. The filter thus has the 
effect of limiting the noise, in some sense, to a band of 
frequencies. Thus we define the effective {or equivalent) 
noise bandwidth of H by 



'Actually R»(0) is the mean-square voltage of the process y(t), so 
the unit of power we are speaking of here is volts’; if a 1-ohm re- 
sistance is assumed, the power is in watts. However, for the most 
part, we shall be dealing with power ratios, such as sign"’ '.o- noise, 
so the units are u n i m po rt ant, as long as they are consistent 
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In these terms, the output noise power is 

N = N 0 W„ |H(M|*. X - (2-25) 

Viewed graphically (Fig. 2-3), XV h is the spectral width of 
an ideal band-pass filter whose maximum response is the 
same as that of H(s) and whose output power (in the 
presence of white noise) is the same as that of H(s). 



Some authors prefer to consider a single-sided behavior 
of filters. They argue that negative frequencies are not 
“really” observable, and hence that all the concepts involv- 
ing frequency should be defined accordingly, to involve 
only positive frequency. This is not efficient mathemati- 
cally, for one immediately rules out such powerful analysis 
tools as the ordinary Fourier transform. Hence, these 
authors are led to a double-sided mathematical analysis 
and a single-sided interpretation. 

For example, the filter response in Fig. 2-4 has noise 
bandwidth W = 2B, wherea„ most engineers would agree 
that the band of frequencies passed is only B wide. 
Working with a single-sided spectral system, as in Fig. 2-5, 
one folds the negative frequency response onto the 
positive. The result is the single-sided spectrum, which 
we give a different notation: 

1 2 S„(/w) 

C ”«->= 0 -<o. <2 ' 26) 





The single-sided equivalent bandwidth is also given a 
new symbol 



and in this new notation, the noise power spectral 
density is 


G n »(fa) — 



S 2N 0 for a> > 0 

0 « <0 . 


(2-28) 


The resulting noise output is the same as (2-25) but ap- 
pears written in single-sided parameters as 


N=N,B„ • 


There is really no practical advantage to be gained from 
a single-sided treatment once the engineer realizes that, 
since physical filters have both positive and negative fre- 
quency responses, it is not possible to measure only the 
positive side or onlv the negative side of a spectrum by 
using linear filters — to the contrary, it is only possible to 
measure the composite effects of both positive and nega- 
tive frequencies with a filter. On the other hand, one can 
certainly build a device that finds the autocorrelation 
function of x(t) and then takes its Fourier transform. 
With such a device, negative frequency spectra ate sepa- 
rately observable; however, the result is a minror image 
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of the positive frequency spectrum and hence yields no 
new information. We shall use both concepts more or 
less interchangeably in this work, and notation will make 
it clear which is meant in any particular case: S(/w) is a 
double-sided spectral density, and W is a double-sided 
bandwidth, while G(/w) refers to a single-sided spectral 
density, and B is a single-sided bandwidth. 


2-D. Sinusoidal Filter Inputs 

Let a sinusoidal carrier having power P be inserted into 
the filter H. We represent this input as the stationary 
process 


x(t) = (2P)*cosM + 0) (2-29) 

where 0 is a uniformly distributed random variable, and 
P is a constant (namely, the carrier power). Then we have 


Rxx(r) — P COS <i>0T 

(2-30) 

S**(/a>) = irP [8(<u (i> 0 ) + 8(w + <oo)] • 

The factor of ir is present in S„ because we must satisfy 

i r* m 

Rxx(O) = P = I S xx {jm)du). (2-31) 

By our formula (2-17), the output spectrum is 

SM = irP [8(o>— <i>o) + S(<i> + a»o)] (2-32) 

and the output signal power (call it S) is 

S = R m ( 0) = P |H(/a>o)p. (2-33) 

If white noise is also present at the input, so that 
x(t) = (2 P)* cos (u> 0 t + $) + n(t), then the output con- 
sists of a signal component with power S and noise with 
power V: 

RJfi) = S + N. (2-34) 

The resulting output signal-to-noise-power ratio (SNR) is 

(2-35) 


„ _ S _ P |H(K)| 
P * N N 0 W h ]k(jJ)\ 


2 

max 


If the carrier frequency is placed at the filter's maximum 
gain point, the SNR is maximized, and 


PlTdnu) N 0 W h ' ( 2 ‘ 36 ) 

Throughout this work, we shall consistently use the 
symbol p to denote signal-to-noise ratios; p„ then denotes 
the signal-to-noise ratio of the waveform y\t). 


2-E. Fiducial Bandwidth 


The ratio W w |H(/<o)| ^ ax /\H(jo> 0 )\ 2 appearing in (2-35) 
occurs very frequently in the theory of phase-locked loops, 
as we shall see in later chapters. It is therefore very 
convenient to give this quantity its own special notation 


w H = 2b H = 


1 I"*™ 

W H \ H(M|U _ vL |H(M| 

wm? w 


2 dc 


(2-37) 


This quantity is much like the noise bandwidth of H, and, 
in fact, reduces to W H when the maximum filter response 
occurs at o> 0 . It is the same form of definition as that given 
to W H , except that it is referenced to an arbitrary fre- 
quency <i>o rather than to the filter's maximum response. 
For this reason, wc shall refer to it as the fiducial band- 
width of H(s). This fiducial bandwidth is then the spectral 
width of an ideal band-pass filter whose response at <■> — w 0 
is the same as that of H(s) and whose output power (with 
a white noise input) is the same as H(s). 


2 -F. Band ’Pass Mixers 

Now consider what happens when a sinusoidal signal 
in white noise is heterodyned (a nonlinear operation) into 
a different passband by a device such as that shown in 
Fig. 2-6. Let the premixing filter H(s) have bandwidth 
W H = 2B h , as shown in the Figure. The filter input, 
denoted by x(t), is 

x(t) = (2 P)* sin (<i > 0 t + 0,) + n t (t), 

where n,(f) is white, with spectral density N 0 , and P is the 
power in the sinusoid. The premixer filter output y(t) is 
then of the form 

y(t) = (2P) V4 |H(/io 0 )|sin (<o„t + 0 2 ) + n 2 (t). (2-38) 

Here 0 2 = -f arg H(/a>o), the new noise n 2 (t) is band- 
limited by H(s) to the filter width W H , and S„ 2 „, (/«>) 
= N 0 \H(jo>)\*. 

We shall assume now that y(t) is multiplied by a unit- 
power 3 sinusoid, 2* sin (u k t + 0 3 ). Thus, the mixer out- 
put 4 is 

v(t) = P*\H(M \ COS [(<a 0 — w*)t + (0 2 ~ 0j)] 

— P vk |H(;«) 0 )| COS [(«1) 0 +O)*)t + (0j + 0a)3 + «a(t). 
(2-39) 

This assumption can be relaxed by including a multiplier gain K m 
in the output v(t). 

This supposes an ideal multiplier. Physical realization of a multi- 
plier may leave a carrier term and harmonics. 


( 
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•Ji COS ( u t +Q j) 



The bandpass filter F(s) following the mixer is centered 
at a frequency <«i, which can either be a> 0 + m or <u 0 — 

For convenience, say it is the latter, <ux = o» 0 — “h , with 
to 0 ^ m- One normally chooses B F ^ B H , since there is 
nothing to be gained with B F > B„. The new noise term 
n,(t) is also split into two bands, and unless the band- 
widths of F(s) and H(s) are chosen properly, overlapping 
of the noise spectra will allow unnecessary noise in the 
output of F(s). It is easy to see from Fig. 2-7 that this 
can be avoided by choosing 5 

2 Bf + Bh < 4 fk 

Bf < 2f h or W F < 4 f k • (2-40) 

When this is satisfied, the output z(t) from F(s) is 

z(t) = P» | H(Jmo) | • | F(M | cos ( Wl t + 9<) + 

(2-41) 

This may need to be modified by a factor of Vi if carrier is present 
in the multiplier output. 


with 04 = 02 - 03 + arg F(j» i). The output noise spectral 
density is, for positive values of «>, 

S» 4 „» = |n 0 | H [#•+«*)] |* • | F(M |* . 

(2-42) 

We have tacitly assumed that F(s) is chosen in accordance 
with (2-40) to suppress completely the image term that 
would appear in S„ 4 „ 4 (i<*>) centered at <i> 0 +<d*. 



Fig. 2-7. Choosing filter bandwidth* to avoid the imago 
noise problem 
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The equivalent noise bandwidth of the mixer is 


W* = 


J_ f” 

"Jo 

\H[i(u,+m)] | 2 *| 

F(;<o) | 2 d<x) 


«[«- + ■*)] F(b) |ma* 


— W F . 
(243) 


The output signal-to-noise ratio of the mixer is then 
given by 


Pz N ( ,W U | W[;(o,+ WA )] F(j<j>) | ^ ax • 


(2-44) 


resulting output occurs at baseband (i.e., o, = 0), and 
the post-multiplication filter F(s) is now a low-pass, rather 
than a band-pass, device. No longer are there four sep- 
aiate bands of frequencies in the multiplier output v(t) 
as shown in Fig. 2-6. Instead, the two low-frequency 
portions merge into one at zero frequency. Because of 
this, theie is an apparent decrease, by a factor of 2, in 
the detector bandwidth, for we note that a band of fre- ! 
quencies W H = 2 B„ at carrier frequency (see Fig. 2-6) 
is heterodyned to a width WVt = Vi W H — B H at base- 
band. At the same time, the merged noise spectra add 
together. 


This attains its maximum value if H and F are built so 
that 

I H(M I = I H(M I ' 

and 


The sinusoidal signal components of the input are also 
heterodyned to zero freouency, and depending on the 
phase of the mixing frequency, these can add or subtract 
in varying degrees. With wj = 0 in Eq. (241), we have 

z(t) = P* | H(ju> 0 ) | F(0) cos 0 4 4- n 4 (t) (2-46) 

where 0 4 = d 2 — 0 3 . Because of the merging of the het- 
erodyned bands, the spectrum of n 4 (t) becomes 


in which case 


Pi (max) 


p 

n»<wi~iwi 

2 

max 

N 0 W U 

JH[/( U +,»)]FW | 

2 

max 


P P 

N 0 Wj, ~ N 0 W F ■ 


(2-45) 


S« 4 , 4 </.) - ^{|f/[/(«-^)] | 2 + |ff(^,+- # )|*}|F(^)|* 

- No | H[/(«-*> 0 ] | 2 1 F(/co) | 2 . (247) 

The latter approximation is valid when F(s) is much nar- 
rower in bandwidth than is H(s). The output signal com- 
ponent thus has a zero-frequency component with density 


This is one of the important features of a simple 
product mixer: the output signal-to-noise ratio is the 
same as that of a simple sinusoid passing through a linear 
filter having the same bandwidth as the mixer. Thus, 
even though mixing is a nonlinear operadon on the in- 
coming waveform, there is no degradation in the signal- 
to-noise ratio at the output. Also, there would be no 
nonlinear distortion if a bandlimited signal were put into 
the mixer, rather than a sine-wave. The heterodyne merely 
produces a translation in the frequency domain and may 
be treated as a linear device otherwise. Care must be 
exercised to avoid image noise, and this can be done by 
choosing B F < 2 f h (balanced mixer) or B F < f h (non- 
balanced mixer). 


S«,.,g> (/«>) = P | ff(K) 1 2 F 2 (0) E( cos 2 0M- 

(248) 

The term E( cos 2 0 4 ) above is referred to as the coherence 
factor. When 0 2 and 0 3 are completely unrelated, then 0 4 
takes on any value between 0 and 2 * with equal likeli- 
hood, so the coherence factor in this case is 

E( cos%) = . (249) 

But when fi s and 0 3 are fully coherent, i.e., when 
0 t ~ 0 3 with probability one, we have 

E( cos s 0 4 ) = 1 . (2-50) 


2-G. Amplitude and Phase Detectors 

When the heterodyne frequency is exactly equal to 
f 0 , the mixer becomes an amplitude detector or phase 
detector, depending upon the relative phases of the in- 
coming wave and the mixing signal. In either case, f he 


Thus, there is a 3-db advantage in signal-to-noise ratios 
to be gained by properly phased synchronous detection. 
In u ay event, the output signal-to-noise ratio is 

- — E(cos 2 0 4 ) 

max 

(2-51) 


0 p I H(M | 2 F 2 (0) 

Pz ~ No W D | H[/(o,-W)]F[;M] 


io 
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where the detector bandwidth W D is given by the same 
formula (2-43) as was the mixer bandwidth W*. By 
choosing H and F to have maximum response at /„ and 
0, respectively, the maximum SNR is achieved 

P P 

Pn max) NnW[} '^ 1 * * 4 ) ~ 7V 0 W/,' ^( C0 * ”^ 4 ) * 

(2-52) 



■o 


S ¥¥ {j^)-ZkTR 


o 



The last two equations resemble those of the mixer, 
viz., (2-44) and (2-45), except for the coherence factor. 
Also, one should realize that the passband of frequencies 
around the carrier is 2W» wide, whereas the detector 
bandwidth is only W hence, care must be taken not to 
confuse the carrier passband with the detector bandwidth. 

For a synchronous phase detector, of course, $ 2 and d. t 
are purposely made to differ by an angle of about ir/2, 
so very little signal power is present in the output. The 
signal-to-noise ratio of a phase detector output is thus 
not a very’ meaningful quantity. 

All of the expressions for signal-to-noise ratios above 
could have been more simply expressed in terms of 
fiducial bandwidths. The treatment here, however, illus- 
trates the desirability of locating the signal power at the 
maximum response point of the filter. 

2-H. Noise 

Many sorts of electric signals are called noise. In the 
early days of radio, noise was familiar as the crash and 
crackle of static. Later, there was the rasp of ignition 
noise and the hiss of thermal and shot noises generated 
in radio circuits themselves. In the end, many engineers 
have come to regard any interfering signal of a more 
or less unpredictable nature as noise. 

The study of noise began with the consideration of 
certain physical sources of noise and the sorts of noise 
that they generate. At first, only very simple properties 
of the noise signals so generated were understood and 
described. As the art has progressed, a mathematical 
theory of noise has grown up. This theory is a part of 
the general field of statistics, and it deals with signals 
that have an unpredictable, statistical, random element. 

1. Thermal Noise 

Perhaps the most fundamental noise is Johnson noise, 

the noise from a resistor. The engineering fact is that a 
resistor of resistance R acts as a noise generator with an 

equivalent circuit as in Fig. 2-8. In an ordinary resistor, 


Fig. 2-8. Thevenin and Norton equivalent circuits of 
noisy resistors 

it is a summation of the effects of the very short current 
pulses of many electrons as they travel between collisions, 
each pulse individually having a wide spectrum. In this 
case, the noise is a manifestation of the Brownian move- 
ment of the electrons in the resistor. In a resistor consist- 
ing of two opposed, close-spaced, hot, electron-emitting 
cathodes, it is a result of the current pulses of randomly 
emitted electrons passing from one cathode to the other. 
In a lossy dielectric, it is the result of random thermal 
excitations of polarizable molecules, forming little fluc- 
tuating dipoles. 

For any particular sort of resistor, it should be possible 
to trace out the source and calculate the magnitude of 
the Johnson noise, and indeed, this approach has been 
used. 

Any energy involved in thermal noise must clearly 
come from the surroundings in the form of heat transfer. 
Thus, derivations of the behavior are destined to involve 
thermodynamical arguments. 

Consider a network containing many resistors. If we 
heat one hotter than the rest, energy tends to flow from 
the hot resistor to the cooler resistors. Johnson noise is 
such energy flowing as electric power. Even when the 
resistors are all at the same temperature, power will 
flow back and forth between them through the connect- 
ing network, always so that, on the average, a resistor 
receives just as much power as it sends out. 

Statistical mechanics tells us how much energy must, 
on the average, be associated with each degree of free- 
dom of a system when the system is in thermal equilib- 
rium. In an electrical network of inductors, capacitors, 
and resistors, the number of degrees of freedom is the 
number of inductors plus the number of capacitors. 
(Inductors in series and capacitors in parallel are treated 
as a single inductance or capacitance, because in setting 
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up a signal on the network, we are free to specify arbi- 
trary initial currents in all the inductors and arbitrary 
voltages across all capacitors.) 

Classical statistical mechanics says that in a system 
(in our case, such an electrical network) that is in equilib- 
rium (all at the same temperature T) there is, on the 
average, an energy VikT joules associated with each 
degree of freedom: 

k ~ 1.380 X 10" 23 joule/degree Kelvin. 

According to quantum mechanics, the energy is less 
than this at high frequencies; Nyquist and others have 
used the quantum-mechanical expression to get the cor- 
rect result. However, even up to many thousands of 
megacycles, the classical expression is accurate. The only 
change that quantum mechanics makes is to say that 
the mean energy per degree of freedom is 

iA4‘ T+ - (2 - 53 > 

where h is Planck’s constant 



E(i z )-kT/L E(v z )-kT/C / 

Fig. 2-9. Total Johnson noise current- and voltage- 
squared relations in simple circuits 

Let E(i 2 ) be the total mean-square noise current in 
the inductor. We can write 

±LE(i*) = ±kT. (2-54) 

On the left we have the average power in the inductance. 

On the right we have the average value this must have 
according to statistical mechanics. Accordingly, 

E(i 2 ) = (2-55) 


h = 6.625 X 10~ 34 joule seconds 

and / is the frequency. This reduces to the classical 
result when f« kT/h. At 300°K 

^ = 6.25X 10 12 cps 


This must be true regardless of the value of R. If R is 
low, we have a narrow-band circuit; if R is high, we have 
a broad-band circuit. 

In a similar way, in the case of the capacitance C and 
the conductance G in shunt, we easily find that 


Hence the theory holds, to engineering accuracy, up to 
60 Gc. Even at very high frequencies, one may define 
an equivalent noise temperature T y as 


T n 


hf/k 

ghf/kT _ J 


= T + . . . 


In this way the energy per degree of freedom is always 
kTjf/2 (anyway, T* is usually inferred from noise power 
measurements rather than calculated from 7'). Johnson 
noise thus often serves as a reference for the noisiness 
of radio receivers and amplifiers, and the effective noise 
temperature provides a useful way of specifying the 
noise output of any source. 


Let us consider two simple circuits as particular exam- 
ples, to see how things work out. In these circuits an 
inductance L is in series with a resistance R at tempera- 
ture T, and a capacitance C is in shunt with a conduct- 
ance G at a temperature T, as shown in Fig. 2-9. 


icE(v’)=±kT 


(2-56) 


If the conductance is small, we have a narrow-band 
circuit with high low-frequency noise components. If 
the conductance is large, we have less low-frequency 
noise but more bandwidth. 


The relations described above of course apply to 
capacitors and inductors not merely in the simple circuits 
we have considered, but to capacitors and inductors any- 
where in all circuits, no matter how complicated they 
may be. In any case, we see that the noise voltage- or 
current-squared is proportional to the temperature T. 

Among the circuits to which (2-55) and (2-56) apply 
is the simple RLC series resonant circuit of Fig. 2-10. 
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c 





The value of the integral can be computed by substitut- 
ing (w/o>o = e* to reduce it to a standard form in integral 
tables; it is equal to 



Qi)d(w/ 0>o) 



ir . 


(2-61) 


We thus have a white spectral density of the equivalent 
noise voltage in the circuit: For every a> 0 


S V ,(M = 2 kTR . (2-62) 

This is commonly referred to as Nyquist’s Law. 


Fig. 2-10. Tuned circuit for investigating the spectral 
density of v n (t) 


The spectral density of the equivalent shunt noise 
current in a resistor can be computed similarly. 


This circuit is characterized by a resonant frequency f 0 , 
a Qo, and an impedance Z(/«), related by 

1 

“° (Lcr 



(2-57) 


1 z 1 2 - «• + (- L - tb)’= [> + <x(t~ t)]. 

We can regard this circuit as one excited by the Thevenin 
equivalent Johnson noise voltage v n (t). If we make the 
Q 0 very high, so that the bandwidth of the circuit is very 
narrow, we can find out something about the spectral 
density of this voltage. Clearly, we must have the voltage 
and current spectra related by the equation 


Sa(jot) 


|Z(/ W )K 


(2-58) 


The integral of Su(/<o) yields E(i 2 ), given by (2-55). 
Thus 



At very high values of Q 0 , supposing S v n vjju>) is fairly 
constant about <u 0 , we can bring S VnV Jju>) outside the 
integral: 


kT 

L 


S " n vJM 

2irLR 


f* foo 

I r 


Qgd(u)/ tl>o) 


+ <?; 


(— — 2 
\<0<l 0 ) ) 


(2-60) 


What happens if we connect two resistances in series 
or two conductances in parallel? In a given frequency 
range, the voltages or currents produced by different 
resistances are uncorrelated; they have random phases, 
and the mean square of the sum of the separate voltages 
or currents is equal to the sum of the mean-square volt- 
ages or currents of the separate resistors. 

As we have noted, for a complex impedance the series 
noise voltage generator at any frequency can be calcu- 
lated from the resistive component fl of the impedance, 
and the shunt noise current generator from the conduc- 
tive component G of the admittance. We can calculate 
the thermal noise spectral density for any network, 
simply by associating with each resistance a series volt- 
age generator according to the relation above. 

We may ask, what is the thermal noise power N 
available from a resistor? We will draw off the maximum 
power if we supply a matched load of the same resist- 
ance. Thus, the available noise power in the bandwidth 
Wn can be obtained by calculating the noise power 
flowing into a resistance R from a source with an internal 
resistance R and an open-circuit voltage spectral density 
given by (2-62). This power is 

N =yW, (watts). (2-63) 

This result is usually stated in the literature in single- 
sided notation: 


N = kTB„ (watts) . (2-64) 
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The conclusion to be reached from this discussion is 
that when an amplifier is matched to its source, the 
amplifier input noise power spectral density S„ B (/o>) is 

S„„(M = No = ~fi (volts 2 ) 


or 


G nn (/a)) 


N + = kTR for / ^ 0 
0 f< 0. 


(2-65) 


fluctuation in velocity is also responsible for what is 
commonly called the modified or reduced shot noise in 
space-charge-limited flow of electrons from a cathode. 


A very simple theory of noise in space-charge-limited 
diodes and triodes (at frequencies low enough so that 
transit time is not important) predicts that the noise can 
be represented by an impressed noise current fit) in the 
plate circuit with spectral density 


One last point concerning thermal noise: By considera- 
tion of the entropy of a system in thermal equilibrium, 
it can be shown that individual electrons have a Gaussian 
velocity distribution of zero mean in each direction. The 
composite effect is a linear combination of Gaussian 
processes, and hence the total noise voltages in a circuit 
are independent white Gaussian processes. 


S u (/o>) = (0.644)2fcT ( g . (2-66) 

Here T c is the temperature of the cathode and g is the 
conductance of the diode or the transconductance of the 
triode. More elaborate theories lead to a factor that, in 
various circumstances, may be a little greater or a little 
less than 0.644. 


2. Shot Noise and Other Noises in Electron Tubes 

Electricity is not a smooth fluid; it comes in little 
pellets, that is, electrons. The flow of electrons in a 
vacuum tube is accompanied by a noise of the same 
nature as the patter of rain on a roof. Schottky, who 
first investigated this phenomenon, called it the 
Schroteffekt (from shot); it is now usually called simply 
shot noise. Like Johnson noise, shot noise has a fairly 
flat spectrum. This is really what we should expect of 
a random collection of very short pulses (each of which 
has a wide spectrum). However, at frequencies for which 
the transit time from cathode to plate is comparable to 
the period, the noise induced in the plate circuit is a 
function of frequency. 

Random processes other than the random emission of 
electrons can also give rise to noise. For example, noise 
can be introduced into an initially noiseless electron 
flow if the electrons randomly hit or miss the wires of 
a grid, with a certain average intercepts- v» of current. 
Such noise is called partition noise or interception noise. 
If a small fraction only of the current is intercepted, 
the added noise is roughly equnl to shot noise for the 
intercepted current. 

At high frequencies and long transit times the excita- 
tion of a circuit may depend on the velocity of the 
entering electrons. In such a case, the random variation 
of velocity of emission from one electron to another, 
associated with the Maxwellian velocity distribution of 
electrons leaving a cathode, can give rise to noise. This 


At moderate frequencies, noises in triodes and diodes 
agree fairly well with (2-66), being perhaps a little higher. 
At low audio frequencies and below, flicker noise appears. 
This typically, but not always, has a 1/f spectrum (dis- 
cussed below). Flicker noise is very variable from tube 
to tube. It has been ascribed to fluctuations in the work 
function of the cathode surface. 

3. Noise With a 1/f Spectrum 

Johnson noise is in a sense inherently white noise, in 
that the fundamental relation between the noise source 
(the resistance or conductance) and the amount of noise 
per unit bandwidth is independent of frequency. Shot 
noise is, in the same sense, inherently white noise too, 
although it can give rise to different spectra in circuits 
with different transfer admittances or in tubes of different 
transit times. In a close-spaced diode formed of opposed 
cathodes at the same temperature and with no average 
current, shot noise and Johnson noise are two names for 
the same thing. 

Some important sorts of noise are generated only in 
nonequilibrium systems, for instance, in systems in which 
dc current flows. Among these are contact noise, such as 
is produced in a carbon microphone, and the noise pro- 
duced in carbon resistors and in silicon and germanium 
diodes and transistors. Both contact noise and transistor 
noise have a spechum that varies nearly as 1/f over a 
large frequency range, though it may be constant at high 
frequencies and at very low frequencies. One could 
obtain a 1/f spectrum down to any given frequency by a 
proper mixture of pulses of various lengths. 
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Transistor noise has been attributed to the trapping 
of the holes or electrons (carriers) that form the current 
flow. The trapping and subsequent release of a charge 
carrier is equivalent to a rectangular pulse in the current. 
The effect may be strengthened by the charge of the 
trapped carrier modulating the flow of other charges. By 
assuming a particular distribution of trapping times, or 
pulse lengths, a 1/f spectrum can be obtained. The matter 


of noise in semiconductors is by no means thoroughly 
understood, and somewhat different mechanisms have 
been suggested. 

Actually, the power spectrum cannot vary as 1/f right 
down to f = 0, for this would imply an infinite noise 
power. Measurements do show a 1/f spectrum down to 
frequencies as low as 10 "* cps, however. 
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CHAPTER 3 

FORMULATION OF THE LOOP EQUATION AND BEHAVIOR 
IN THE ABSENCE OF NOISE 


In this Chapter, we shall develop the basic equation 
governing the phase-locked loop. Assuming that noise is 
absent, we shall develop the acquisition and steady- 
state character of the loop. This essentially follows the 
work of Viterbi. 


The term 0(f) appearing here is the loop estimate of 0(f), 
and K, is the rms output of the VCO. 

The result of this multiplication, if perfect, would be 


3- A. The Basie Integro-Differential Equation 

The essentials of a phase-locked device are a multi- 
plier, a loop filter, and a voltage-controlled- oscillator 
(VCO), as shown in Fig. 3-1. Most of the more elaborate 
systems using double-heterodyne techniques, IF limiters, 
and acquisition aids reduce to this basic model, insofar 
as analysis of the behavior is concerned. The input is 
assumed to be a sinusoid of the form 

x(t) = A( 2)'* sin [«o„f + 0(f)] (3-1) 

in which the quantities A, <u 0 , and 0(f) are 
A = rms voltage amplitude of x(t) 

<t>o = frequency of the VCO when its input is shorted 
6{t) = the input signal phase process. 

Normally, 0(f) consists of an information-bearing term 
tf/(t) due to modulation and a term d(t) due either to 
doppler shift in the received signal, to drift in the VCO 
center frequency, or to some bias voltage appearing at 
the VCO input. 

The incoming signal is multiplied by the VCO output 
v(t) = K, (2)* cos [to„t + 6(f)] . (3-2) 


t/,(f) = AK i (sin [0(f) — 0(f)] + sin [2to 0 f + 0(f) + £(f)]}. 

(3-3) 

However, generally speaking, the multiplication is ac- 
complished by a device unable to respond to the double- 
frequency term. Also, the multiplying device has some 
gain K m , and hence the actual output of the “phase- 
detector” is 

y(t) = AK t K m sin <f>(f) . (3-4) 

The quantity 4>(t) = 0(f) — 0(f) is called the true phase 
error. 

The multiplier output is then fed into the loop filter 
F(s) and emerges as z(f), which, in turn, supplies the 
input to the VCO. 

In many actual implementations, the waveform y(t) is 
monitored (perhaps multiplied by a gain factor) and is 
called the “dynamic phase-error.” Similarly, z(f) is often 
referred to as the “static phase-error.” These are func- 
tions of sin <f>, rather than the true phase error <f>, and thus 
tend to generate some confusion. Little reference to these 
terms will appear in the material that follows. 


*(/) -Ajz »in(to 0 /+ 0<n) y(t) -AK x K m tin +(t) 




•'(f)**', Sz eo« + 0(f)) 



Fig. 3-1. The basic phate-locktd loop 
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The VCO output frequency is, as its name implies, a 
linear fu* 1 ‘on of its input. 

*>VCo(t) == «»o “t" Kv CO Z(f) "i" Kvco ® • (3*5) 

Consequently, the phase estimate 6(t) developed by the 
loop, being the integral of the VCO frequency, can be 
written, omitting e for the present, 

6(t) = *Z2>z(t) (3-0) 

where p = d/dt is the Heaviside operator. Substitution 
for z(f) yields 

${t) - AK,K m K V co ® sin . (3-7) 

P 

At this point it becomes convenient to define 

K = KiK m K\co (3-8) 

as the open-loop gain ® of the loop, and to substitute 
§= 0 — <f>. This produces the fundamental equation that 
specifies the behavior of the loop in the absence of noise, 

0(t) = + AK ^ sin 4>(t) . (3-9) 

When F(p) is a constant, (3-9) is a first-order integro- 
differential equation. Hence the configuration of Fig. 3-1 
is called a first-order phase-locked loop. Similarly, when 
F(s) has n finite poles, the system equation is an (n + l)th- 
order one, and the device in Fig. 3-1 is said to be an 
(n+1 )th-order loop. 


and, in such a case, the loop error is related to the input 
phase by a simple, linearized version of (3-9), 

m = P - A * F - ( ^ 4>(t) . (3-10) 

Generally, the input process can be separated into an 
information-bearing part i/-(t) and a phase-offset term d(t): 


p + AKF(p) + p + AKF(p) ' 

(3-11) 

We see that there are two kinds of error present in the 
loop. That part of the error due to the modulation <f>(t) 
is commonly called phase distortion. The remaining error 
is that produced, for example, by a doppler shift on the 
incoming carrier. To be effective, the loop must be de- 
signed to track whatever variations d(t) may have, so 
the filter F(s) must be properly chosen. Error arising 
from d(t) is called the tracking error, or transient 
distortion. 


The steady-state tracking error can be determined by 
the final-value theorem of Laplace transform theory. 
Denoting D(s) = X[d(t)], the Laplace transform of d(t), 
this theorem reads 


4>f — 


,. s 2 D(s) 

To s + AKF(s) • 


(3-12) 


An example will elucidate the behavior: 


3-B. Tracking When the Loop Error is Small 

The phase-locked loop would be of little use if it were 
not possible to use it to reconstruct the input phase 
process with some degree of fidelity — that is, unless we 
are able to “lock” the loop in the first place. The term 
“lock” is somewhat subjective at this point, but it may 
generally be thought of as a condition in which <f>(t) never 
varies outside an interval of size 2». The mean-phase 
value is called the lock-in point; such lock-in points are 
located 2* radians apart. Whenever the phase error goes 
through 2w radians, we say that the loop has skipped a 
cycle. If the loop is capable of reducing the phase error 
to a small enough value, say | <f> | < w/6, we can approxi- 
mate 


sin <t> <f> 


'Often the open-loop gain K is defined to contain the dc gain of 
F{«). However, F(0) may not always be finite (e.g., if F(«) has 
an integration), so we do not include it here. In cases where it is a 
finite value, we shall be careful to design F(0) = 1. 


First-order loop: Suppose F(s) = 1, and d(t) = $„ + 
n„# + Vi A„ t 2 , where <t>„ is the initial phese offset, n„ the 
initial frequency offset, and A„ the doppler rate. 


D(») = v + f + -F ■ 


This produces a steady-state tracking error given by 


4>„ = lim 

«->0 


n„ *f A (]/ s 

s + AK 


(3-13) 


The error is not bounded if A„ ^ 0; but when A 0 = 0, 





(3-14) 


In such a case, 4>„ can be made small by adjustment 
of the loop gain. N’ote that the loop will track any initial 
phase offset 6. with no steady-state error in the absence 
of a frequency offset (i.e., n 0 = 0). 
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From this example, one can see that if d(t) = 

(the nth derivative of d(t) is A», so D(s) = A and 
if F(s) has l poles at the origin, as 

m = * °> «°> ** 0 f 3 - 15 * 


then the steady-state tracking error would be 
A - T im A n p(^'-"-‘ 

s Ul p(s) + AKq(s) 


(3-10) 


This is not finite if n > / + 1. When n = l + 1 the error 
is finite and has the value 


<f>n 


_ A, p(0) 

AK q(0) 1 


(3-17) 


To investigate the behavior for larger values of n„, let 
us rewrite the loop equation, differentiating (3-9), 

n 0 — 4> + AK sin ^ . (3-20) 

We shall denote the frequency error by = fl(t), and 
plot n versus according to (3-20). The result is shown 
in Fig. 3-2. Whenever fi is positive, <f> tends to increase, 
and whenever n is negative, 4> tends decrease. Note 
that if | fi 0 1 < AK, there are regular points at which 
n = 0. Starting at <f>( 0) = ntr, where n is an even integer, 
the system tends to move along the sinusoidal trajectory 
of Fig. 3-2 until it reaches the <£- axis at 

= n* r + Sin-‘(n„/AK) . (3-21) 


When n < /, there is no steady-state error, <j>„ — 0, and 
perfect tracking ultimately results. 

The conclusion here is the following: In order for a 
phase-locked loop to track an nth-degree phase function, 
F(s) must have at least (n — 1) poles at the origin. 

3-C. Acquiring Lock in the First-Order Loop 

Let us concentrate now on the way a first-order loop 
(F(s) = 1) behaves in a tracking mode only. We consider 
the case in which (modulation) is absent from the 
input phase function $(t); only the d(t) term due to input 
doppler, VCO drift, etc., is assumed to be present. From 
the discussion in the previous Section, it is clear that the 
first-order loop- can track a function of the type (constant 
frequency offset) 

6(t) = 9 0 + n 0 t (3-18) 

whenever n 0 is small enough that 4>(t) — 0; that is, when 


This is a stable point; (l cannot become negative because 
<f> would then tend to decrease and return the system to 
the ?-axis. If n were an odd integer, the system would go 
through a larger part of the sinusoidal trajectory until it 
reached a stable point at 




= (n + !)*■ + Sin -1 • 


(3-22) 


If n 0 > AK, however, the trajectory never crosses the 
and phase lock is never achieved. The maximum 
pull-in range of a first-order loop is thus 

Om« = AK . (3-23) 


Whenever ft 0 < AK, the loop ultimately tracks the in- 
coming function 6(t) with no frequency error, but with a 
constant lagging phase error, given by (3-21) and (3-22). 
The loop never skips a cycle. 

The multiples of * here are usually emitted, and the 
steady-state phase error is written merely 
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that is, the steady-state output of the VCO is a replica of 
the input, lagging by <f>„: 

Vt»(t) — 2 '* COS [(«,„ + ft,,)t " 4 * On • 

(3-25) 


The pull-in time may also be determined from (3-20) as 
follows: we can write 


dt_ _ 1 

d<f> ft 0 ~ AK sin <t> ' 


(3-26) 


This can be integrated to give the pull-in time to any 
particular value of 4>, but, since the denominator above 
vanishes at the lock-in point, an infinite time is required 
before ft = 0. However, if we agree that the loop is in 
lock whenever | <f> — 4>t 0 ^ | < Stock, the corresponding 
time is finite. The maximum time required to achieve 
lock, designated the acquisition time, will be that 
time required for the system to pass from — * + S| 0 ,. k 
— Sin -1 (ft JAK) to Sin" 1 (n ,,/AK) — 8 tork . When the value 
of 8i nrk is small, the answer r duces to approximately 


t»cq 


AK cos 8 lock 


(3-27) 


l ! -.nee, if a loop is designed so that 4>„ = 5 deg when ft 0 
is 200* rad/sec, and 8 tr ck is taken to be 5 deg, then the 
required pull-in time is approximately 1 msec. 

The trajectory described in Fig. 3-2 is called a phase- 
plane diagram of the loop behavior. It is particularly 
useful for analyzing the lock-in characteristics of low- 
order loops. 


The steady-state phase error <£, s is the same as for the 
first-order loop, 

= sin" (-^jr) , (3-30) 

and, clearly, the loop never locks if ft,, > AK. Even for 
values of ft,, less than AK, pull-in possibly may not occur, 
even though a stable point exists in the phase plane This 
is a consequence of the fact that the trajectory that 
takes tj> to <£„„ is not the simple sinusoid of Section 
3-C, but is rather the solution to Eq. (3-29). 


To simplify things a bit, we can substitute for 4> the 
value 


.. _ rfft dft (U (/ft 

dt d 4> * dt ” <4 


(3-31) 


(note that 4 ,s proportional to the slope of the phase- 
plane trajectory). This produces the trajectory equation 



d ft 
d<t> 


+ 1-1 AKt 2 cos ft + AK sin . 

(3-32) 


A solution to (3-32) depends not only on ft,,, but also 
on ft(0) and the rnitial conditions of the VCO. For 
certain values of ft,,, the loop will lock regardless of what 
values fi(0) and <f>„ take. See, for example, Fit . 3-3. In 
this Figure, <f> has been limited to the region (--*,*) 
by folding all the trajectories onto this region. When ft 
is positive, 4> increases, and when ft is negative, <f> de- 
creases. Hence, motion as a function of time is from left 
to right in the upper half plane, and from right to left 
in the lower. Starting a trajectory, say in the upper half 
plane, one follows it to the right until 4> = », skipping 
back »o <f> = — * at the same value of ft encountered at 
<f> — *. This continues until the lock-in point is reached. 


3-D. Acquiring Lock in the Second-Order Loop 
With Passive Loop Filter 

The system of interest in this Section is one in which 
the loop filter takes the form 

= (3 - 28) 

The discussion in Section 3-B indicates that a loop with 
this filter can follow a phase function of the type 
d(t) — 0,, + ftj (i.e., a constant frequency offset), but 
that no higher-order terms may be considered. The equa- 
tion governing the behavior, from (3-9), is 

fto = tj + (1 + AKr t cos 4 ) i + AK sin ^ . (3-29) 


For values of ft u larger than those in Fig. 3-3, lock-in 
may occur for some initial conditions of the VCO but 
not for others, as illustrated in Fig. 3-4. This latter Figure 
shows that there is a limit cycle toward which all higher 
trajectories converge, as well as some of those from be 
low. Even if ft(0) is negative, lock-in is not assured, since 
only those trajectories that ultimately pass through a 
strip determined by the asymptotes of the saddic- points 
can converge to the lock-in point. 

By integrating (3-32), Viterbi has derived a necessary 
condition on ft v : If lock hi occurs for all initial conditions 
of the VCO, then ft 0 is bounded by 

1 0.1c 2 f (~j (l + KAK„)1 (M3) 


1 
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Fig. 3-4. Lock-in bohavior of a second-order loop with imperfect integrator, FIs) = (1 +r 2 s)/(1 + r,s), for n JAK 
— 0.9 and AKt\/t x = 2. The upper bound of Eq. (3-33) Is QJAK <0.693. Note that lock-in 
occurs only when the trajectory happens to pass through the “slot." Otherwise, 
the trajectory enters the periodic frequency lag region shown. 
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Experimental evidence seems to indicate that this is a 
fairly sharp bound; that is, for n„ very close to the bound 
above, lock always ultimately occurs, although as the 
bound is approached, longer and longer times are re- 
quired to pass through the region where the limit cycle 
would occur should Q 0 reach the bound. 

Acquisition time for this loop is difficult to compute, 
but as an approximation, Viterbi has computed the time 
required for frequency lock (i.e., that time required be- 
fore no more cycles are skipped); when t, is large, it is 



The additional time to acquire phase lock is probably 
about the same as f acq of the first-order loop, given in 
(3-27). 


3-E. Tuning the VCO 

In the first four Sections of this Chapter, we have 
neglected the effect of the VCO tuning voltage e shown 
in Fig. 3-1. Including it in the VCO equation (3-5), we 
can write 


<ovc()(0 — <u 0 + KvcoC "f~ Kvco z(t) . (3-35) 

The effect of e is that it determines the VCO frequency 
when z(t) is absent. The loop phase estimate is likewise 
changed to 

9(t) — ((D 0 + Kv coe)t + ^ r0 z(t ) . (3-36) 

The basic integro-differential equation governing the 
loop is almost the same as (3-9), except for a term K vco et 
subtracted from 6(t). The result can be written 

6(t) — Kvcoe = n(f) + AKF(p) sin </>(<) . (3-37) 

The conclusion is now evident: All the answers we have 
obtained about lock-in, tracking, etc., are the same as 
before, except that we replace 6 by 6 — K vc 0 e. Alterna- 
tively, we can redefine u> 0 as <a op = <o 0 + K vco e. 


For example, -he lock-in points now occur at 





0 « Kvc o8 


AK 


)■ 


(3-38) 


By choosing e properly, the steady-state error can be 
made zero. This value of e is clearly 


e = 


n„ 

Fvco 


(3-39) 


In the first-order loop, the acquisition time with the 
foregoing value of e becomes 



(3-40) 


3-F. Locking the Second-Order Loop with 

Perfect Integrator / 

When t i » Tv, the system whose loop filter is given by 
(3-28) behaves much the same as the one with 

m = . (3-41) 

T lO 

except that the latter is able to track any constant fre- 
quency offset with zero steady-state error. However, 
there is a steady-state error when a doppler-rate A„ is 
present: 

= Sin_1 (^) (3-42) 

and just as there was a limit on 12 0 in Section 3-D for 
lock to occur for all initial conditions of the VCO in the 
passive-filter loop, there is a limit on A 0 in the perfect- 
integrator loop. 


Clearly, from (3-42), lock cannot occur at all if 
A 0 > AK/ t,. Figure 3-5 is a typical phase-plane diagram 
of the way such a loop behaves with a doppler-rate 
input. Almost all trajectories with fi(0) < 0 eventually 
lead to the lock-in point, while most of those with 
fl( 0) > 0 diverge upward, never reaching lock. The Fig- 
ure is drawn with the largest value of A 0 (determined 
experimentally) such that all trajectories with n(0) < 0 
ultimately terminate at the lock-in point, while most of 
the trajectories with fi(0) > 0 fail to lock. This value of 
A 0 is a function of the parameter AKr'- /n = r; 




(3-43) 


Figure 3-6 illustrates the way k(r) varies as a function 
of r. 


When the VCO is being swept as in Section 3-E to 
lock onto a constant frequency offset n,„ the rate by 
which e is changed should not produce a A„ exceeding 
(3-43), and, as Fig. 3-5 shows for positive values of A 0 , 
the sweep should begin with a negative frequency error 
(a positive frequency error if A„ is negative). That is, 


Kvco^l 


de 

dt 



(3-44) 


In the passive-integrator loop, (3-43) and (3-44) hold 
approximately for short periods of time (see Section 
6-C3). 
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* 

Fig. 3-5. Phase-plane trajectory of a second-order loop with perfect integrator to a doppler-rate 

input A 0 for AK t\/t x — 2, and A 0 = AK/ 2r t 
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Fig. 3-6. Normalized maximum doppler rate, Mr) = A„ rJAK, for which 
lock is guaranteed in the absence of noise, as function of the 
loop parameter r = AK r\h\ 
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CHAPTER 4 

BEHAVIOR OF PHASE-LOCKED LOOPS WITH STOCHASTIC INPUTS 


The consideration given to phase-locked devices in 
the preceding Chapter excluded both the possibility that 
the incoming phase process 8(t) was a random, informa- 
tion-bearing process and the possibility that the input 
wave was a noisy one. We now wish to include these 
cases for consideration. Some comments are in order 
concerning the effect of noise on the pull-in and tracking 
behavior under such conditions. A small amount of 
phase jitter due to noise or 6(t) will not affect the phase- 
plane trajectories significantly; it will result only in a 
small amount of hash being superimposed on the tra- 
jectories. Such a small amount of jitter on the trajector- 
ies will thus not appreciably change the number of cycles 
required to achieve lock. On the other hand, if the input 
noise is significant, the disturbance of the trajectory may 
easily be such that a drastically different number of 
cycles will be required before lock is achieved. The 
number of cycles required to achieve lock depends on 
the particular noise waveform observed at the input 
and is random. If there is a great amount of phase jitter, 
the response of the loop may be so erratic that lock never 
seems to occur. There are, as yet, no analytic results on 
the pull-in time when stochastic inputs are observed. 


4-A. Development of a Mathematical Model 
and a Basic Loop Equation 

We start as we did in Chapter 3, but now the input 
consists of a sinusoidal signal with power P = A 2 im- 
mersed in a noise ni(f): 

x(t) — A 2 1 * sin [(o 0 f + 0(f)] + Mi(f) . (4-1) 

The density of Ui(f), denoted S n . n .(/«>), will be arbitrary 
for a moment, although ue shall restrict it later. 

Following the same procedure as in Chapter 3, we 
derive the basic integro-differential equation, similar to 
(3-9): 


- AK m sin , (1) + + K H‘ + "•MJ - 

v ' p r p p 

(4-2) 


Again, <f>(t) — 0(f) — d(t), K = K x K m K v co, e is & tuning 
bias applied to the VCO, and the term n v (t) is the internal 


loop-noise referred to the VCO input. The noise n(f) 
appearing in (4-2) is 

n(f) = n f (t) 2’* cos [o» 0 f -I- 0(f)] . (4-3) 


Except in special cases, the spectral density S nn (j<o) of 
n(f) is not merely a heterodyned version of because 
of the correlation between n,(f) and the phase estimate 

d(t). 


As we indicated in Section 3-E, the effect of e on the 
loop is that it changes the no-error VCO frequency from 
(i> 0 to w 0 p = wn + KvcoC. We can thus merely redefine u>« 
to include this offset and proceed with 

*> = sin m + ™ „(,) + J^L „„(,) 

= sin [0(f) -0(f)] 

+ «(f) + n„(t ) . (4-4) 

This last equation gives us the basis of an exact 
mathematical model of the phase-locked loop. As shown 
in Fig. 4-1, there are three inputs to this model of the 
loop: the first, 0(f), enters, is differenced with 0(f), then 
passes through a nonlinear, zero-memory amplifier whose 
characteristic is y out = A sin x ln . The input system noise 
n(f), added to the output of the sine-amplifier, becomes 
the input to a filter with transfer function KF(s)/s, whose 
output is then $(t); in this way, 0(f) satisfies (4-4). 

Solution of (4-4) is really made no easier by the model 
in Fig. 4-1; however, we are not generally interested in 
the solution to (4-4) for a specific sample function n(f) 
anyway. What we are interested in is some measure of 
the loop’s average performance, such as rms phase error, 
mean rate of cycles skipped, etc. From the Figure, how- 
ever, one is led to some natural approximations that aid 
in loop analysis, and these are developed in the ensuing 
chapters. 
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INPUT 

nit) 



*+ n ¥ ( t) 


Fig. 4*1. Exact mathematical equivalent of the phase-locked 
loop. Sources of external, as well as internal, noises 
are shown; VCO tuning voltacje e is also indicated. 


4-B. Discussion of Mathematical Models 

As a first approximation, the output phase-jitter of a 
phase-locked loop can be obtained by replacing the 
sinusoidal nonlinearity of our mathematical model (Fig. 
4-1) by a linear or polynomial approximation. However, 
the steady-state-error probability density p(<f>) is a periodic 
function because the loop is sensitive only to errors 
modulo 2tt; consequently, the steady-state phase variance 
in every phase-locked loop is infinite. This is borne out 
by the fact that, in actuality, the loop skips cycles at a 
certain mean rate, executing a random, nonstationary 
motion between lock-in points, much like a discrete ran- 
dom walk. The linear and polynomial approximations 
mentioned above do not exhibit this periodic lock-in 
behavior, and hence they not only give no information 
concerning loss of lock, but even come up with a finite- 
variance, stationary phase-process* a case that never 
actually occurs, as we have said. 

But while the steady-state phase process itself does 
not possess a finite variance, it does when phase angles 
are reduced modulo 2ir. On this basis, Tikhonov and 
Viterbi were able to provide exact results for the first- 
order system. 

Considering the exact equivalent model, there is cer- 
tainly no loss in generality by the explicit inclusion of 
the “mod 2*” nonlinear function in the loop as in Fig. 
4-2. However, there are several differences one can inject 
into the mathematical results: Although the steady-state 
^-process is nonstationary, the <f>(mod 2w)-process usually 
is stationary (assuming 9(t) is stationary); there then 
exists an equivalent stationary process (call it 4>(f) ) with 


u lite variance, which, if reduced modulo 2ir, would result 
*. a a random process with the same probability density as 
And so long as our sole interest lies in <f >( mod 2ir), 
and noi in c\- auvior of <f>, we may assume that <1> is at 
work in the loop, and that stationarity prevails. In fact, 
the linear and polynomial approximations are direct 
attempts at calculating the behavior of <h. 

4-C. Spectrum of the Loop Noise 

The noise we deal with in our analysis is not that 
occurring at the input, but 

n(t) = n,(t) 2 ,/a cos [w„t + #(<)] (4-5) 

in which $(t) is derived, in part, from nj(f). The non- 
independence of 6 and n, means that S„„(/a») is not gen- 
erally the convolution of S n . ni (j(o) with the spectrum of 
2 H cos + 6). However, it has been reasoned that 
whenever n f (t) has a much greater bandwidth than §(t), 
and whenever S„ jnj (/«>) is symmetric 7 about <o 0 , then the 
convolution formula is approximately correct. This gives 

S„»(/fl)) = S„j B 4 [/(<0 — <i) 0 )] + *£ S n .„.[/(w + oi 0 )]- 

(4-6) 

That is, the loop noise has the same spectral shape as ni(f) 
does, except that it is heterodyned down to baseband. The 
spectral shape of Hj(f) usually comes from allowing white 
noise with density N a = kT/2 to be passed through a 


’This condition is stronger than it needs to be. See Davenport and 
Root, p. 162 . 
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predetection filter H($), whose brndwidth is, say, W/,. 
That is, «i(t) has a spectrum of thi form 

Sn i » 1 (M = (4-7) 

and, in this case, the loop noise n(t) has as its spectral 
density 

S-(M =-|n.{|H [/(«-«.)] |* + |H[/(-+-o)]|*}. 

4-8) 


In the zone about zero frequency, the region of interest, 
these two halve-, add: 


S„„( 0) = N a \H(1*o)\*. (4-9) 

(Again we note the change in bandwidth by a factor of 
one-half due to heterodyning n,(t) to zero frequency, as 
in 2-6.) 


n(t) 



Fig. 4-2. Equivalent exact mathematical model of phase-locked 
loop, with explicit reduction of (mod 2w). The VCO has been 
replaced by ff V ro/s (tuning bias and VCO noise omitted). 


REFERENCES FOR CHAPTER 4 

1. Jaffe, R. M., and Rechtin, E., "Design and Performance of Phase-Lock Circuits 
Capable of Near-Optimum Performance Over a Wide Range of Input Signal and 
Noise Levels," IRE Transactions on Information Theory, Vol. IT-1, pp. 66-76, 
March 1955. 

2. Davenport, W. B., Jr., and Root, W. L., Random Signals and Noise, McGraw-Hill 
Book Co., Inc., New York, N. Y., 1958. 


r - 


27 




! 


JPL TECHNICAL REPORT NO. 32-819 


CHAPTER 5 1 

THE LINEARIZED ANALYSIS OF PHASE-LOCKED SYSTEMS 


Most of the important behavior of a phase-locked de- 
vice can be predicted by replacing the loop phase- 
detector by an ideal phase-differencer. The distinction 
between the two is, of course, that there is a sinusoidal 
nonlinear characteristic for the phase detector, whereas 
the ideal phase-differencer is a linear error-detecting 
device. To justify this linearized analysis, one must as- 
sume that the loop error is very small, which usually 
means that a high signal-to-noise ratio prevails, and that 
modulation on the carrier is not excessive. The latter of 
these restrictions can be removed by subcarrier modula- 
tion. Foi a discussion of such techniques, see Volume II 
of this work. 


As with any linear filter, L(s) has some effective noise 
bandwidth W L , which can be found by using (2-24), 

~2n / I I 2 ^ (5-4) I 

W /v = 2B l = — • ?- V ’ 

where we use L 2 to denote | L(/<o) | 2 ia ^ 

It is important to remember that \V L is the bandwidth 
of a transfer function at baseband and it is not the same 
as the width of the carrier passband produced by the 
loop. As we saw in Section 2~F, the carrier passband has 
width 2 W L . 


5- A. Behavior of a Linear Loop 

If the level of n(t) is sufficiently low, and if the loop 
is designed properly, the phase error <f>(t) should be very 
small. In such a case, the approximation 


A slight rearrangement of (5-2) relates the phase error 
to the input processes: 

4,(t) = [l-L(p)]6(t)-L(p)2@- (5-5) 


sin 4> = <f > , (5-1) 

when inserted into (4-4), yieids a linear equation relating 
the loop input and output 

omitting VCO noise. 

The output phase is thus the result of a linear filter 
acting upon the input phase process 6(t) immersed in a 
normalized version of the loop noise, with the normalizing 
factor in this case being equal to the rms signal ampli- 
tude A. An equivalent circuit of the linearized loop 
appears in Fig. 5-1. The overall loop transfer-function 
(call it L{s)) is related to the loop-filter F(s) by the rela- 
tions 


(omitting for the present the VCO noise term n„(f), to 
be considered later). The first term of (5-5) is an error 
due to the incoming phase function, and, hence, it repre- 
sents distortion due to filtering. There are usually two 
effects comprising 6(t); one is the “doppler” phase shift 
d(t), and the other is the information process f(t). The 
distinction between d(t) and \p(t) is that d(t) is a non- 
stationary, more or less deterministic phase variation, 
whereas is a stationary, zero-mean random process. 
The last term in (5-5) is phase error due to the presence 
of noise at the input of the loop. The mean-square phase 
error, which we shall denote by S 2 , is thus composed of 
three terms: 

2 2 = p?{t) + S 2 + a 2 . (5-6) 

The first term represents the transient distortion due to 
d(t), the second is modulation distortion, and the third 
is the mean-square phase-noise. The latter two, being 
stationary, can be computed by using (2-19): 



1 f* M 

8 !== M„ I 1 - WI’SwWA* 

(5-7) 

, _ N 0 W l L 2 _ N.BlL* 

° A 2 ~ A 2 ‘ 


i 

I 
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| 


% 


The expression for r is a very interesting one; it states 
that the mean-square phase-noise is precisely equal to 
the noise-to-signal power ratio in the loop effective noise 
bandwidth times the squared maximum loop response. 
This must not be given the wrong interpretation. It 
should be realized that N„ is the noise density at carrier 
frequency, that A-’ = P is the incoming signal power, 
also at carrier frequency, and yet W/, is the equivalent 
noise bandwidth at baseband, computed by (5-4)! There- 
fore a 2 is not the noise-to-signal power ratio in the pass- 
band about the carrier frequency, which has width 2 W L , 
as discussed earlier. 


It is convenient to include in d(t) the initial phase offset 
0„. There is usually a complete lack of knowledge about 
6 n , so it is assumed to be a uniformly distributed random 
variable over ( — w, ir), thus with variance w 2 /3. The re- 
mainder of d(t) (call it di(t)) is a time-varying phase 
function (usually due to doppler shift) whose form is 
fairly well known. By Parseval’s theorem, the total tran- 
sient distortion can l»e computed in the frequency domain 
by using the relation 

t = 1 ~ L(h) i 2 i d(m i *] do> • 

(5-10) 


Equation (5-8) can also be written in terms of fiducial 
bandwidth Wl (or b L ) as 


NqWl W 

A 2 A 2 2 


(5-8) 


and our previous statement can be amended to read: the 
mear-square phase-noise is precisely equal to the noise- 
to-signal ratio in the fiducial loop bandwidth. 

The total transient distortion is defined by the equation 

" S(t) dt . (5-9) 



Note the resemblance this bears to $ 2 in (5-7). The term 
D(s) is the Laplace transform of d(t). Equation (5-10) 
results when one substitutes E [D(s)D(— «)] for S^(s) 
in (5-7). 

We shall have occasion to investigate these quantities 
more fully in Section 5-C, to determine the best loop 
configuration. 

5-B. Calculation of Loop Bandwidth 

The loop bandwidth formula (5-4) is rather easy to 
apply to the loops with simple filters. For future refer- 
ence, we shall tabulate a few of these. 
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1. First-Order Loop 


The corresponding value of L 2 is given by the equation 


The first-order loop is one with no loop filter, F(s) — 1. 
Hence (5-4) gives 

to, = Wl = \ AK or = = (5-11) 

This relation enables one to express the loop transfer 
function as 


L(s) = 


2W L 

s + 2W L • 


(5-12) 


U 


2 ( i+ ? + 2 r )X 1_ ^)" + ( r4 ^)"' 2<r+1) 

T 2 

~2(l+2r)* + r 8 -2r-2 ' (5 ‘ 16) 

This last function (Fig. 5-3) is monotone decreasing in 
r, approaching unity as p-*cc, and becoming infinite as 
r-* 0. 


2. Second-Order Loop, Passive Integrator 

The second-order loop in most widespread use is one 
in which the loop filter takes the form (see Fig. 5-2) 


m = 


1 + T 
1 + T\$ 


L(s) = 


1 + TtS 

1 + (r, + l/AK)s + (r ,/AK)s’ 


(5-13) 


Then, from (5-4), the equivalent loop noise bandwidth 
can be computed, 


W, = ^ - 

2r-j (1 + Ti/fT\) L*’ 

_ (r + l)[r- - 2r - 2 + 2 (1 + 2r) y ‘] 

2^P 


(5-17) 


(again we have normalized F(0) = 1, as agreed to in 
Chapter 3). In the usual case, F(s) is designed with 
ti » Ti, to make F(s) appear to have a pole at the origin. 

\ 

It is convenient to define a quantity r as the following 
ratio: 


the latter approximation being valid when rr, » r 2 . The 
fiducial bandwidth (2-37) is a somewhat simpler quantity: 


Wl = 


r + 1 

2r* (1 + Tj/fri) 


r + 1 

2 . 

2 r 2 


(5-18) 


AKr| 

*1 


(5-14) 


The value of u producing maximum loop response, 
U = | L(fu) |* M , can be found by differentiation; it is 



[(1 + 2r)*— 1]*. 


(5-15) 


O 1 w. 


o 



F(t) * 


I + r 2 » 
l + r,J 


WHERE r 2 * R X C 


r, ~-(R y < C 


Fig. 5-2. Passive integrator loop filter 



Fig. 5-3. Variation of maximum loop response P, noise 
bandwidth Wl, fiducial bandwidth w L , and frequency 
at maximum loop response, for second-order phase- 
locked loop, as a function of the parameter 

r = AKtJ/t, 
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The loop transfer function is then easily expressed in 
fiducial bandwidth: 



From (5-13), one may compute the system damping be- 
havior. For r < 4, the roots of the system are under- 
damped, with damping factor f and loop natural 
frequency /3 given by 

C = i-(14 r -j- 






(5-20) 


For r = 4, critical damping (( = 1) occurs, and for r > 4, 
L(s) has two real negative poles. 


3. Second-Order Loop, Perfect Integrator 

Whenever an operational amplifier is permissible in 
the loop, the loop filter takes the form 

F(s) = ^p- S . (5-21) 

TiS 

We have kept the notation here the same as that in the 
previous exar*ple for later convenience. This loop filter 
gives the same response as that given in (5-18), (5-19), 
and (5-20) when the r 3 /rr, term is dropped. This indi- 
cates that whenever rr, » r s , the loop with an imperfect 
integrating filter (5-13) performs very much the same as 
that with the perfect integrating filter (5-21). 


J-C. Optimization of Loop Paramoters 

According to the preceding sections, we can perform 
a linear analysis to determine the behavior of a phase- 
locked device whenever there is sufficient justification to 
warrant the assumption that sin ^ But whenever 
noise is present, there will be arbitrarily large phase 
deviations if only we wait long enough; because of this, 
the assumption above cannot be strictly valid at every 
instant of time. 

It is natural, then, to design the loop transfer function 
to minimize the effect of noise, that is, to choose K and 
F (s) to minimize the mean-square ’"op error. In this v-a y. 


not only does the loop operate in the most hnear manner, 
but also 6(t 1 is most faithfully reproduced by 8(t). As can 
be seen from (5-6), the mean-square phase error is not 
independent of time; to minimize 2^ thus would require 
a time-varying loop filter. 

Rather than derive an optimum time-varying filter, 
Jaffee and Rechtin modified the criterion somewhat, bas- 
ing their optimization on minimizing a somewhat differ- 
ent quantity. They define the total phase error, denoted 
by 2* , as 

2 * = A* « ; + 8 s + <T» (5-22) 

in which is a Lagrange multiplier, a design parameter 
related to the bandwidth of the loop. 


Minimization of 2£ by choice of the loop transfer func- 
tion L($) specifies both AK and F(s). But the question 
now is, “How does one choose L(s) to minmize 2J?" 
Jaffe and Rechtin recognized that A-cj, -I- 8 2 could be 
computed by replacing S^(s) in (5-7) by A •E[D(s)D(—s)] 
+ Sw (s). This led them to the conclusion that is 
minimized whenever L(s) is chosen in accordance with 
the Weiner optimization technique, which yields the 
Yavits-Jackson formula 


L„ p , (s) — 1 


N„'*/A 

TW 


(5-23) 


S(s) = A- E [D(s) D(- a)] + S++ (s) + N 0 /A* 


Some explanation of the bracket terminology here is in 
order. The bracket [ ]‘ refers to a type of “square-root" 
factoring of the enclosed function, retaining in [ ]* the 
left-hand-plane poles and zeros of the enclosed function 
only; singularities on the imaginary axis are equally 
divided between [ ]* and its mirror image [ ]\ 

The mechanics involved in (5-23) will be made clearer 
in the next Chapter. 


5-D. Tho Effects of VCO Noiso 

One of the limiting factors governing the design of 
narrow-band phase-locked devices is the phase noise 
Inherent in the output of the VCO. This noise appears as 
a random fluctuation, drift etc., and is sometimes called 
“oscillator bistability." More often, however, the term 
“instability" is used in context to mean “frequency insta- 
bility,” whereas the quantity of concern to us here is the 
phase stability of an oscillator. 


I 


I 
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To represent this effect, let n,(f) be that noise voltage 
applied at the input of a perfect noiseless VCO which 
produces the same noise as appears in the actual VCO 
output. That is, is the VCO phase noise referred to 
its input. The closed-loop behavior of the loop is now 
somewhat modified by the addition of this internal noise: 

6(t) = L(p) [*(*) + -^j + [ i-W ] Kvco n v (t) . 

(5-24) 

There are many factors that contribute to the spectral 
makeup of n r (t). The two most significant terms that 
appear are (1) thermal (Johnson) noise generated in the 
resistances of the oscillator, and (2) noise with a 1/f 
spectrum, as discussed in Section 2-G, associated with 
the transistors, varactor diodes, carbon iesistors, etc. 
Thus, the spectrum of Ky C on v (t) can be approximated by 

&VCO (fa) — N 0v + 2nN lv / | w | . (5-25) 

As we have previously indicated, the 1/f law cannot 
ext 2 nd all the way down to zero frequency, so the equa- 
tion above is really valid only when j w j is greater than 
some small value c. But if [1 - L(s)]/$ has a zero at the 
origin, this c need not be known. 



Fig. 5-4. Factors governing relative contributions of VCO 
noise to output phase noise 


1. Optimization of w L and r 

Depending on the values of N 0 /A 2 , N or , and N lv , there 
are optimum choices for both w L and r. 

As shown in Fig. 5-4, there is quite a broad range 
of r for which g (r) is fairly constant and is nearly equal 
to its minimal value of 1.5491 at r = 5.22. A quite useful 
range to use this approximation for g(r) is from about 
r = 3 to r = 10. Outside this region, g(r) is increasing, 
drastically when r < 3, and more slowly for r > 10. 


The amount of phase error in the closed-loop output 
due to VCO noise can then be found by integration: 

+00 

-- 1 7 [No, + a r N u A] (L . 

(5-26) 


On the other hand, the coefficient of N or , i.e., (r -h l)/4r, 
is a monotone decreasing for all r, asymptotically ap- 
proaching the Vcdue 0.25 as r becomes infinite. If r is 
larger than 5 or 6, for most practical purposes, we can 
use the value 0.25 with little fear of producing any sig- 
nificant error. 



At this point, for our treatment, we shall assume that 
L(s) is the P **ssive-inte ^ator loop of (5-13), and again 
set r = AKt;/- l T^e tr .al phase error due to input and 
VCO nois* in ^s of the fiducial bandwidth w L , is 


„• = ^ g(r)-^f . (5-27) 

The function g(r) is shown in Fig. 5-4; it is given by 

J*. + *)• + W ~ lilfor r > 4 

4[r* (r -- 4)] (r — 2 — \/r(r — 4) 


,1 J 25 

*W = 'U6 


(r 

\2[r°(4 


-i ) 2 r* 
-»•)]*[ 2 


tan - 


=±\ 


Vr(4 - r) J 


for r = 4 

tor r < 4. 

(5-28) 


As a conclusion, then, we see that the best value of r 
to minimize a- lies to the rigbl of 5.22 and probably is 
less than 10, if any N„. is present at all. In fact, the dif- 
ference between <r; nn and tr foi any r between 5 and 10 
is almost inconsequential insofar as tr is concerned. The 
same statement cannot be made for r < 4 (an under- 
damped loop). Thus we may take a value of r = 7 as 
being practically as good as r 0P t, regardless of the other 
parameters. (This reasoning alleviates the necessity for 
differentiating a 2 to find the exact value of r opt .) 

We can now use this value of r, differentiate <r 2 with 
respect to w L , equate to zero, and solve for w L to find its 
best value. Straightforwardly, Wl is the solution to the 
equation 

(5-29) 
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2. An Example 

Just to see how a typical design should be made, let us 
assume that we have parameters 

~ =6X10* 

No 

N or = 0 (5-30) 

N„ = 0.08 . 

The first of these is typical for an Earth-spacecraft link 
at a distance of about 3 million kilometers, while the 
values of JV,„. and N u . have been approximated from 
plots of a' 1 under very high signal-to-noise conditions. 

The optimum value of tv L is then 
= 26cps 


and the optimum loop transfer function is 
r ( s \ - 1 + 0-154» 

opt w 1 + 0.154s + (3.38 X 10~ 3 ) s- ' 

These parameters also produce a phase deviation of 
a = 2.5 X 10' 2 rad = 1.43 deg rms . 

3. Conclusion 

It should perhaps be mentioned at this point that de- 
signs of this type are most valid for oscillator “clean-up” 
loops and spacecraft carrier tracking loops, or in situa- 
tions where the tone to be tracked is spectrally a very 
pure one. Design of the ground receiver tracking loop 
would probably not use this analysis, since its loop must 
track phase deviations imposed on the carrier by the 
noise in the spacecraft system. There are other factors 
that must be considered, such as frequency acquisition 
interval, lock-in time, and doppler tracking rate. What 
has been presented here should be taken merely as a 
guide as to what the ideal bandwidth is from a minimal 
noise point of view. 
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CHAPTER 6 

OPTIMIZED DESIGN OF TRACKING FILTERS (LINEAR ANALYSIS) 


A phase-locked receiver can be used either to track 
the cairier component of the incoming signal, or to 
demodulate its information, or both simultaneously. In 
space communications, the amount of doppler shift on 
the Cc.rrier is always measured, because it provides infor- 
mation about the spacecraft velocity. A carrier tracking 
loop is thus always present in a space-communications 
receiver. 


6- A. Tracking Loop Design 

Suppose it is our aim to provide the best filter to track 
a given doppler-phase polynomial d(t) of degree N- 1, 
assuming for the present that modulation of the carrier 
is absent. The form of D(s) is then 

DM = 7 + f+-=^ <e-i) 

in which the degree of Q(s) is less than N. The filter 
specified by (5-23) is 

Lnpt(s) = 1 — j= TTaTav * 

( 6 - 2 ) 

In the examples worked by Jaffee and Rechtin, the 0 O 
term was always set equal to zero. This corresponds to 
the case in which the loop is initially tracking with no 
phase error. 


6-B. Optimum Filter for Random Phase Offset 

The simplest example of loop optimization occurs for 
d(t) = 0 O > a uniformly distributed random phase-offset. 
To find Lovi(s), wt insert N = 1 into (6-2): The denomi- 
nator in (6-2) is 


| n , A 2 A‘ J it 2 ~ | + _ , n\A f£} 

s * + N„ 3 J “ s + (3N„)^ ‘ (6 ‘ 3) 


Thus, the optimum filter is given by 




XAir 


S + 


J 

XA7T 

(!W.)' 1 J 


(6-4) 


The loop bandwidth, specified by (5-11), is 


Wj 4 — W L = 


\Att 

2(3N,r ‘ 



We can now eliminate the Lagrange multiplier to give 
the optimum loop design equations in terms of the loop 
bandwidth, viz., 


L apt(s) ~ s + 2 { V/ 

F(s) = 1 (6-6) 


K = 2 WJA . 


6-C. Optimum Filter for Frequency and 
Random Phase Offset 

The next example of interest is the optimization of a 
loop with frequency offset Q 0 : 

d(t) = + n 0 1 (6-7) 


where again is random. With N — 2 and @(s) — 
O 0 + fto s inserted into (6-2), the denominator of L op ,(s) 
becomes 



A 2 AV , A 2 A 2 Q 2 "| + 

3 N„ S N 0 


s 2 + 


/A 2 AV , 2>.An n Y A AAn„ 

V 3N„ ~N? ) N? ' 


( 6 - 8 ) 


The loop natural frequency is tliu.- 


- = A AQn 

H \7 Vfe 


(6-9) 


in which case the optimum filtering function is given by 


s 

2/8 2 + 

^ 3n 2 

V4 

+ |8 2 

s 2 + 

V- .l s±' 
L ^ 3n 2 

]?+*’ 


(6-10) 
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Comparing (6-10) with (5-13) and (5-14), we recognize 
the optimum parameters: 


T-i 



7T T 

«» [3r(r-2)]“ 



i2i (r - 2) 


r 


2 + 


3Q, 1 ; 


(6-11) 


io h = 


( t + l)n„ 

2nr 


[3r(r — 2)] ' 


We have thus translated the effect of the Lagrange mul- 
tiplier implicit in fi to the design parameter r. It is 
important here to note that r must always exceed 2 in 
optimized loop design. This comes about because of 
the assumption that was a uniformly distributed ran- 
dom variable. 

The Jaffee and Rechtin example using 0„ = 0 produces 
a result somewhat different from that above: 


t 2 = 2 ,/ 7 0 



S Jaffee-Rechtin. (6-12) 

1. Choice of Parameters 

In most cases, fin is not known prior to locking the 
loop (or else the loop would have been pretuned so as 
not to encounter any frequency offset in the first place). 
Hence fi n must be treated as a random variable. An ex- 
amination of the average transient error in (5-9) reveals 
that is minimized whenever E(Q„) = 0 (that is, the 
loop should be initially tuned to the expected incoming 
frequency) and that the fi„ of our previous calculations is 
replaced, in this case, by the rms frequency offset, 

Once suitable values of A/N'f, tv L , and fi 0 have been 
established, /? (and hence r) can be found, and this speci- 
fies what L„pt(s) is to be used. The optimum loop filter 
F op ,(s) is related to L opl (s) by (5-3): 

F„pt (s) = ±±7^ • ( 6 ‘ 13 ) 

TiS 


But rather than synthesize this which would re- 

quire an operational amplifier, it is quite usual to substi- 
tute a passive filter whose characteristic is approximately 
the same as (6-12); for example, whenever » r., 

F «w~fr^ < e - U) 

can be used. This will introduce a steady-state phase 
error (Fig. 6-1) 



which does not exist with the F„,„(s) in (6-13). Once the 
loop is locked, however, the VCO can be retuned to 
eliminate 

As a further consideration, one cannot expect a very 
good lock-in behavior when Q„ is so large that the carrier 
frequency falls outside the initial loop passband. 

On the other hand, it is usually desirable to design a 
tracking loop with its bandwidth much narrower than the 
initial frequency uncertainty region. Such loops are usu- 
ally frequency-swept (by controlling the VCO input 
voltage e ) slowly through the uncertainty interval to 
acquire lock. It is thus reasonable to design the tracking 
loop with parameters to ensure that lock-in proceeds 
optimally whenever the carrier enters the fiducial 8 loop 
passband, i.e., when Q„ = 2,irb L . Upon inserting this con- 
dition into (6-11), the proper value of r can be computed 
numerically: 


r = 2.28245 
£ = 0.755 . 


(6-16) 


This compares very favorably with the results obtained 
by Jaffee and Rechtin (6-12). The optimum loop param- 
eters for r = 2.282 are 


t 2 = 1.643/tCi, 


r a /AK = 


1.180 


(6-17) 


ard the corresponding optimum loop transfer function is 
given by 

r , , _ 1 + (1.643 /w L )s 

opt w 1 + (1.643M> + (1.18/u>z, 2 )s 2 ‘ 

(6-18) 


“Purely for computational convenience here, we have specified the 
fiducial bandwidth hi : rather than the noise bandwidth Bt, as 
defining the edge of the passband. 


39 





JPL TECHNICAL REPORT NO. 32-819. 




(c) FIRST-ORDER LOOP, FREQUENCY STEP 





(d) PASSIVE INTEGRATOR SECOND-ORDER 
LOOP, DOPPLER-RATE INPUT 





Fig. 6-1. Response of first- and second-order loops to various inputs 
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2. Evaluation of Transient Error 

The total transient error is given by the expression 


steady-state error is r — 1; that is, the loop is an un- 
derdamped one with damping factor ( = 1/2. 



(6-19) 


X [1 -L.pt (a)] [1 -L. pt (-a)] ds . 


With the form of loop filter given in (5-21), the transient 
error is 


— 


TiTzQl 

2AKr 


r \ L + 30;T iy / ' 


For the optimum loop, then 

2 * 3 (r ~ 1) 


6n„ [3r (r - 2)T 


( 6 - 20 ) 


( 6 - 21 ) 


With the value r = 2.282, and Q 0 = -kWl, 


( t (opt) 5.37 /%Vi 


(6-22) 


The parameters in the Jaffec-Rectin example (6-12) pro- 
duce almost the same transient behavior: 


f T (JR) 


= 5.4/tc,. . 


(6-23) 


3. The Effects of Doppler Rates in Second-Order Loops 

If we assume that there is a slow doppler rate A„ 
(rad/sec 2 ) superimposed on a simple doppler shift, the 
discussion in Chapter 4 indicates that the passive-filter 
(5-13) second-order loop cannot be expected to maintain 
a small steady-state error. The ultimate phase motion 
due to this doppler rate, A 2 = A„ in the notation of (3-17), 
is linearly expanding in time, as 


_ A n , A 0 (r + l) 2 / t \ 

<^«A(dopp rate) AK * ~ 4 rto 2 \u) • (6-24) 


On the other hand, when the perfect integrating filter 
second-order loop (5-21) is used, there is a finite value of 
the steady-state error, viz., 


Another observation that can be made is the following: 
In using a passive-filter second-order loop to track a 
signal having a doppler rate, it is necessary to compensate 
periodically for the ever-growing steady-state phase 
error. This is conveniently done by manually retuning the 
VCO at least every sec (as discussed in 3-5) to zero 
the phase error due to doppler rate. In this way (6-24) 
never exceeds (6-25). 

4. Comments on the Choice of r in Second-Order 

Tracking Loops 

In Section 5-D, we concluded that a value of r between 
6 and 10 is desirable to minimize the effects of VCO 
noise. In Section 6-C1, we concluded that an underdamped 
loop with r = 2.282 is needed to minimize the total phase 
error. And in Section 6-C 3, an underdamped loop with 
r = 1 provides the best tracking of a signal with nonzero 
doppler rate. 

Normally, in a communications system, the contribution 
of VCO noise is not a factor of utmost concern. On the 
other hand, minimizing VCO noise is of great importance 
in designing frequency-generator “clean-up” loops. Too, 
the doppler tracking capability of a loop is critical in the 
specification of many receivers. Some choice is thus avail- 
able to the design engineer in picking an r to suit his 
particular application. 

5. Case With VCO Noise Included 

Comparing (5-7), (5-10), and (5-26), one can see that 
VCO noise effects can be included in the optimizing 
analysis quite easily when S» t , Ml , (s) is a rational spectrum. 
However, for the case we have considered, the VCO (in- 
put) noise has r l/f term, or more precisely in the s-plane, 
N lv /K 2 co (— s 2 ) l/a , and (5-2.?) and (6-2) are not valid 
when such a branch point appears in one of the spectra. 
A Wiener optimum L(s) certainly exists, but the resulting 
loop filter mav not be a rational function of .v. It is more 
meaningful to limit the discussion here to optimization 
under the constraint that F(s) be rational. 


_ A.Ti 

tpsMdopp rate) \ r* 


A „(r + l) 2 
4rwj t 


(6-25) 


There is naturally a great resemblance here, and at t = r u 
(6-24) and (6-25) are equal. 


In either of the two cases above, it is evident that, with 
other things being equal, the value of r causing the least 


Roweve, we easily recognize that the optimum loop 
for an initial frequency-step offset again must be of 
second order, and that the loop filter must be of the same 
type as that given in (5-21). Then, too, we know that, 
depending on values of the various parameters, the 
optimum value of r lies between 2.282 and about 7. It 
remains only to find specific values of r and to*,. 
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We have computed separately all of the terms appear- 
ing in the total phase error £*, (except distortion due to 
carrier modulation). Hence 

2 2 = A 2 * 2 + <7 2 

_A 2 r 1 T 2 n 2 / , SAK\ N U W L 

2 AKr \ L ^ 3Q 2 ri ) ' A 2 

+ ( r + 4. S( lM± l ' (6-26) 

4 r tv i, tvl 


Applying the philosophy that fi„ be set equal to ntv L 
(so that lock-in proceeds optimally once the offset fre- 
quency enters the loop passband), we reduce this to 

v- = + AV ( r + 1 ) !t 

Zr A 2 16 r 2 tv L 


(r + 1) / A V 
4 rtv L y 3 



gjr)N lv 

wi 


(6-27) 


The value of A, of course, is a parameter left unspeci- 
fied. Previously (in 6-9), we related it to the loop natural 
frequency fi (and through r, to the loop bandwidth, as 
well), which then became the arbitrary parameter in the 
loop. We can also do this in the present case; inserting 
(6-9) into (6-27), we find 


V 2 

“ T 


N i) tVf, / 1 | 1 

A 2 \ 1 + r 


+ 


4r \ 

3 (1 + r) ; * j 


I 


(r + 1 N ar + g(r)N„. (6-28) 

4r tv i , tv’i 


With this particular evaluation of A, there is an optimum 
way to pick r and t v L . However, we can note that none 
of the coefficients, as functions of r, varies drastically over 
the range 2.282 to 7; this means that 2 2 is not extremely 
sensitive to r. Such is not the case with tv L , and the de- 
signer should, in all practicality, seek to determine that 
value of w L which will minimize 2 2 . 
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CHAPTER 7 

THE DESIGN OF SYNCHRONOUS-DETECTOR AGC SYSTEMS 


One of the fundamental circuits built into almost every 
receiver today is some sort of closed-loop regulating 
system that automatically adjusts the receiver gain to 
maintain a constant output level. The most common cir- 
cuit used in broadcast receivers is one in which the IF 
signal is rectified, filtered, and applied to the grids of 
variable-/* tubes in the IF amplifier. Such a circuit (often 
referred to as Automatic Volume Control, or AVC) is 
equally as sensitive to noise voltages in the IF passband 
as it is to signal voltages, and, consequently, affords only 
moderate stabilization of the signal amplitude. 

In phase-locked receivers, it is more usual to use a 
synchronous amplitude detector followed by a very 
narrow-band filter to derive the feedback voltage. The 
output of this filter, chiefly due to the presence of signal 
and relatively insensitive to noise, can be used to control 
the receiver gain very efficiently. 

Such a loop we shall distinguish from the AVC loop 
mentioned above by the name Automatic Gain Control 
or AGC. 

In what is to follow, we shall present a linear feedback 
theory for the design of AGC loops. As a result, one may 


predict performance analytically as a function of loop 
parameters within a certain measure of accuracy. The 
treatment here essentially follows the work of Victor and 
Brockman, with some minor extensions. 

\ 

7-A. The Synchron^us-Amplitude-Detector 
AGC Loop \ 

The block diagram sl\own in Fig. 7-1 shows the prin- 
cipal elements involved in the design of an AGC System. 
Briefly, the input 2 1 * A(t) Vos [<i> 0 f + 0(f)] + n«(t) enters 
a variable-gain amplifier whose gain is a function of a 
feedback control-voltage c\f). The output (amplitude- 
stabilized), appearing as 2* \ A cos [<a<,f + 0(f)] + M0> 
is passed to a phase-locked 'loop to derive a coherent 
reference, against which it is synchronously detected. At 
this point, an external gain-reference voltage ea is com- 
pared with the detector output; \the resulting error volt- 
age, filtered by Y(s), then control^ the receiver gain in a 
way that tends to null any differences between eg and 
the detector output. \ 

There is a definite relation between; AGC voltage c(f) 
and the receiver gain (call it 1/A*\ r c(f)], or merely 
1/A*(f)), for if the receiver output is\held nearly con- 
stant at an rms value of A, it follows \hat c(f) is also 

\ 

\ 


x 0 U) = A(t)JZ cos 



Fig. 7-1. A synchronous-detector AGC loop, using a phase-locked loop to provide a coherent reference. 
(The actual phase-locked loop may be part of the gain-controlled receiver.) 
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related to the rms incoming signal amplitude A(t). Thus, 
c(t) can be monitored by a voltmeter as a calibrated 
measure of input signal strength. 


be the adjusted value of receiver attenuation at a particu- 
lar operating point if c(t) could be set equal to zero. The 
actual receiver attenuation K rec corresponding to K R is 10 


When double-frequency terms of the detector output 
are omitted, the input to Y(s) is 

ZaocW = ~e a + AK d cos <f> (t) + K D ni(t) 

i Mt) v eo _ in\ 4. Kpno(t^ 

- e » + cos + A*(t) ' 

The constant K D refers to the gain of the synchronous 
detector, but otherwise the notation follows that previ- 
ously used. 

Victor and Brockman recognized that if the AGC loop 
is performing well, one may approximate 9 


—e a + ‘fifi K D cos*~e 0 ln ; 
A*(t) l 


A(t) Kp cos <f> ~ 
A*(t) e 0 


e 0 In A(t) - e 0 In +e 0 ln cos <j>(t) . 

(7-2) 


Krec = K H + 


20 '<*(§)• 


Now define the quantities 11 

a(t) = 20 log A(t) 


a*(t) = 20 log — . 

&.D 


(7-4) 


That is, a(t) and a*(t) are the signal strength and adjusted 
receiver attenuation expressed in decibels. Since In x = 
log x/log e. 


c(t) = KcY(p) 


+ (i$i) logOTS * (,)+ ^r] 


_ o*(t) - K r 

K a 


(7-5) 


The error here is less than half the square of the left-hand 
quantity. 


The next step in the analysis is that of choosing a func- 
tion for the variation of receiver gain with AGC voltage. 
Mathematically, it is convenient to express the gain in a 
logarithmic Taylor series: 


epA*(t) _ jq[*, + * a c<o+...]/ 2 '> 

Kd 


20 log £ 


g 0 A*(f) ~ 

K d 


= k r + KAt) + 


(7-3) 


From this equation, we can draw the equivalent circuit 
shown in Fig. 7-2. Note, however, that the loop is linear in 
the logarithms of input and response, rather than these 
quantities themselves. The equivalent loop gain K A oc and 
AGC-loop transfer function C(s) are given by 


Kagc — 


KaPqKc 
20 log e 


C(s) = 


KaooY(s) 

[l + K aqo Y(s)]K a ' 


(7-6) 


Straightforwardly, a(t) and a*(t) are related by the 
linear operator equation 


A characteristic with only the first two terms is nearly 
valid for many voltage-controlled-gain amplifiers, al- 
though Kr and K A depend on the operating point to a 
certain degree. The quantity A*(t)e 0 /K D is an adjusted 
attenuation factor. By retaining only the first two terms 
above, we assume the receiver has an adjusted attenuation 
in decibels which varies linearly with control voltage (over 
a limited range): there is an attenuation of Kr db in the 
absence of control, and an additional K A db/v attenuation 
when c(t) is applied. It should be stressed that Kr would 


'Here we use "In x“ to denote the natural logarithm of x, and 
“log x” to denote the base-ten logarithm of x. 


a*(t) = K A C(p) 

X ^a(t) + g° o 8e ) *° n ‘(*) + 20 lo « 005 *(*)J 

+ [l-K A C(p)]K 8 . (7-7) 


,0 The value of K P . e can be set by a bias added to c(t), if desired. 
However, the equation for K r .« shows this is not necessary, as 
adjustment of eo will do the same thing. 

"Since the input signal power is related to A{f) by ? = AVR », 
where Hi. is the input resistance of the receiver, a(t) differs from 
the value of P in dbw by the term 10 log Ri»: 10 log F = fl(t) 
— 10 log Ru. 
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Fig. 7-2. Equivalent diagram of an AGC loop. Linearized analysis follows the 
assumptions that A*(f) is exponentially related to c(t) and that alM«*a*lf). 
The input is alt) = 20 log Alt). Adjusted loop gain a*(t) = 20 log [A*(t) e o/Ko] 


The AGC voltage is thus related to the input signal level by 
c(t) = C(p) | 'a(t) + Kofii (i t ) 

+ 20 log cos <f>(t) — . (7*8) 

7-B. AGC Stability 

The steady-state mean and variance of a*(t) due to input 
noise are readily found: if A(t) changes slowly with respect 
to the time constants in C(s), 

E [a*(#)] = K A C(0)E[a(t}] -K„C(0)(10 log e) <r 2 
+ [1-K^C(0)]K* 

<7 ' 9) 

+ 2(10 log e)*o^^Jc*(0)K*. 

In this calculation, we have approximated In cos by 
— i <j>*, and we have assumed that is a Gaussian 
process. Both of these assumptions are somewhat in error, 
but behavioral indications are sufficient for our purposes. 
It should be noted that N 0 is the spectral density of the 
stabilized noise n<(t), and not the input noise n 0 (t); for this 
reason, it is referred in (7-8) to A’, the stabilized signal 
level, as N 0 /A J (which is the same as N 00 /A 2 (t) at the re- 
ceiver input, however). As a fair approximation, A = e 0 /K D 
can be used in evaluating the variance of a*(t) above. 
When Kaoc is high, C(0) - 1/Ka . 


The fluctuation of a*(t) depends not only on the input 
signal-to-noise ratios in the phase-locked loop bandwidth 
o>£ and in the AGC loop bandwidth w c , but also upon the 
value of stabilized signal level at receiver output. The 
mean value of a*(t) is less than that of a(t) by about 4.35 db 
at«r 2 = 1. 

7-C. Calibration Equation 

Equation (7-7) shows how attenuation varies with signal 
strength and noise. Moreover, a(t ) is related to c(t) by (7-8). 
Thus, it is possible to relate the steady-state mean value 
of a(t) (call it a), to that of c(t), or merely c, and thus, to 
relate c to the average input signal 12 strength ai 

a = K k + -c(jj) c + ( 101 °g e ) 

= K r .e + 20 log {eo/Ko) + (c|o)) c + (10 log e) a 2 . 

(7-10) 

This calibration curve is well demonstrated in actual 
practice, as attested to by the results shown in Fig. 7-3. 
The receiver whose characteristic is depicted is the S-band 
receiver at Pioneer Site, Goldstone Tra. ’’ing Station, Cali- 
fornia. 

The variance of this estimate of a for a particular 
measured value of c is the same as the variance of o*(1)/Ka 
given in (7-9), because, as shown in (7-8), c(t) and a*(t)/K A 
have the same fluctuations due to noise. 


"If 10 log R i. is subtracted trom the right-hand side of (7-10), the 
calibration curve relates c to signal power in dbw. 
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Fig. 7-3. Measured AGC curve shewing departure from 
linear behavior 


7-0. Dynamic Behavior of AGC Loops 

In calibrating the ACC loop, we assumed that Aft) was 
more or less constant — or was at least changing very slowly 
with respect to time constants in the loop. But this may not 
be a realistic assumption during the early postlaunch 
phases of a spacecraft flight, or even later when the space- 
craft is perhaps tumbling through space. In such cases, 
one must investigate the transient behavior of the loop 
in more detail and, if one can predict what types of 
transients are to be expected, choose parameters to 
optimize performance. 


Fluctuations in signal amplitude can be of two types: 
first, there can be a nonstationary waxing and waning of 
the signal, such as is due to the effects of changes in range; 
and second, there can be superimposed o*\ the first a 
stationary perturbation such as might occur as the space- 
craft antenna seeks to retain its bearing toward Earth, or 
as a result of tumbling. 

The modified Wiener formula (5-23) can be fitted to 
this problem by expressing 


Aft) = 

aft) = 04 ft) + otft) 


where w - have put 


aift) = 20 log A/t) 
a+ft) = 20 log *V(0 


(7-11) 


(7-12) 


and in which we have designated Ayft) and Aift) as the 
stationary and nonstationary parts of Aft), respectively. 
For convenience, a,pft) can be taken to be a zero-mean 
process by lumping the mean value of aft) into a d ft). 
The Wiener optimum filter is 


C.pt(s) 


1 A J E [5d(s)fl<i( '$) J +• Sayayfs) 

“ TsWF * L [*«]- 


pr 

(7-13) 


S'fs) = A \E[S d fs)difs)] + S^*) + (20 log e) 2 (^*) 
+ (10 log e) 2 S^tfs) . 


For convenience, we have approximated e (1 /K D = A in 
the equation above and set & d fs) = X [<**(<)] . 

The bracket [ ]■ is the right-half-plane image of [ ]♦, 
and ( ] pr refers to the physically realizable part of the 
enclosed function. 


The iitter can be computed as XX' 1 [ ] (i.e., the 
Laplace transform of the inverse Fourier transform of 
the quantity). An equivalent method for computing j *,. 
when the enclosed quantity is a rational function is ou- 
tained by expanding the function in a partial fraction 
expansion, but retaining only those terms having left- 
half-plane poles, poles at the origin, and poles at infinity. 


For simple transients, much can be learned about the 
form of C 0 ,t(s) without resorting to solution of (7-13). For 
example, if we are concerned with gain-tracking an input 
of the form aft) = ajft) = A,fVnl, then by a reasoning 
similar to that in Section 3-B, we find (neglecting the term 
due to *'oise) 


E [ fl< (0 - a*(0]«~ 


A ^ O )**- 1 

K AQ cqf0)fn-l)\ 

0 


for n > I 


torn < l. 
(7-14) 


We assumed here that Yft) takes the form 

rw=^. m *o, *omo. (7-is) 


Thus, for stable AGC operation, Yfs) must have at least as 
many poles at the origin as the degree of ajft). 
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As a particular case, when a,i(t) = a, Y(s) may take the 
form 


“ T+^ < 7 ' 16 > 

in which case there will be a constant gain-tracking 
transient error 


E [a„(t) - a*(t)]„ = . (7-17) 


The AGC Ijop bandwidth in both cases is 


w 0 


1+Kaoc 

2t*OC 


k AO( . 

2taoO 


(7-19) 


By making K A oc large enough, a signal-strength change of 
a few db/sec can be tolerated for a short time. 


A better Y(.») to follow a x t would, of course, be of 
the form Y (s) = l/r AO o s, in which case there would be 
a steady-state gain-tracking error 


Bui if 0 ^( 1 ) = a x t, the loop filter (7-16) produces a steady- 
state transient error that grows linearly in time: 





A comparison of (7-18) and (7-20) immediately shows that 
when K A oc is large, the two AL> J loops perform a’ jut the 
same for t < r A ac- 
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CHAPTER 8 

THE DOUBLE-HETERODYNE PHASE-LOCKED RECEIVER 


Communications requirements for deep-space missions 
require receiving systems whose capabilities at first sound 
so extreme as to be unrealizable. The development of the 
double-heterodyne receiver with phase-locked carrier 
reference, however, has made coherent communications 
the standard language of the systems engineer. Indeed, 
except for a few minor refinements, parameter optimiza- 
tions, etc., there are few questions concerning the behavior 
of these receivers that cannot be readily answered. 

A double- i eterodyne receiver, briefly, is merely one in 
which there are two separate intermediate-frequency 
amplification stages. Such receivers can be designed to 
operate at very low signal levels with a great measure of 
stability, precision, and reliability. 


8- A, Basic Configuration of the Receiver 

The receiver shown in Fig. 8-1 combines the advantages 
of intermediate-frequency amplification to produce high 
gains (and there are two such stages) with those of the 
phase-locked loop: namely, coherent communications 
capability at low bandwidths and predictable stability. 
It is not immediately obvious that the concatenation of 
front-end mixing, subsequent IF stages, and ultimate 
phase detection constitutes a phase-locked loop as we 
have seen it in previous chapters. There are practical 
necessities, as well as subtle advantages, to the use of a 
limiter prior to phase detection, and although we have 
nowhere accounted for such a device in our previous 
theory, there is little, if any, deleterious effect due to its 
presence. 

The tracking portion of the receiver operates as follows: 
The input function is 2* A(t) cos [<« 0 f -I- 0(f)] 4- n 0 (t), 
where A(t) is a slowly varying rms signal amplitude and 
n 0 (t) is a wide-band noise with density N 00 . The input is 
mixed with a multiplied-up version of the VCO output, 
so the first IF output process is 

Kn [2* A(t) cos(wjf + 0 — ft) + . 


The heterodyne operation followed by linear gain does 
not alter the signal power/noise-spectral-density ratio, as 
we have discussed in Chapter 2. The first IF frequency is 

, aij 2ir 

2ir too oj/ti 

(We assume u>„ > <o*i as a practical matter, although our 
analysis does not require it.) The signal is now hetero- 
dyned against a locally generated, free-running reference 
of the form cos (&>* 8 (f) — 0,), so the output of the second 
IF amplifier is 

K/ 1 K /2 { 2* A(f) cos [(*>, — o>* 2 )f 4- 0— -0 + 0i] 4- * 

Again, the signal-power/noise-spectral-density ratio is 
preserved, now centered at the second IF frequency 

U ~ 2ir 2ir 

(In general, f 2 may take a negative sign when « > «>i. 
This is a purely mathematical distinction, and need cause 
us no great concern here.) 

Each of the IF amplifiers has a known bandwidth, so 
the overall IF bandwidth or predetection bandwidth 
can be computed as in (2-43); however, the first IF nor- 
mally has much wider bandwidth than does the second, 
so W* is essentially the equivalent noise bandwidth W„ 
of the second IF filter H(s). 

As a matter of theory, either or both of the IF amplifiers 
can be gain-controlled by an AGC voltage, although it is 
practically necessary to control only one of them, say the 
first. Hence, in the terminology of Chapter 7, K n K, 2 
= 1 /A*(t), so the second IF output is, in our previous 
notation, 

x(t ) = 2* A cos (u> 2 1 + 0 — § + Oi) + tii (t) . 

This voltage feeds a band-pass limiter whose saturation 
limits are ± l v. The limiter mean-square output is, of 
course, constant at a value I 2 , of which 2(21 /it) 2 lies in the 
frequency zone about f 2 . A signal component is present 
whose mean-square value is, let us say, 2(2 la/ir) 2 ; the 
factor a 2 is then a signal-power suppression factor, with 
0 < a < 1. The remainder of the limiter output in the IF 
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(7) x 0 {f)zjz 4U>cos[w o / + 0(>)]+/>o</) 

(?) • 'M~'TZ cos [«*, / + B it )] 

(?) *xtit) - Ki 1 1\/2 eos [ {mo - U) A , ) / + + 

( 4 ) K 2 (/) ~ Sz cos (u h2 t-6 ,) 

(?) *(/)= Jz A cos [(w 0 - “*i " «<» 2 ) / + fi- fl + ®|J + /»/(/) 
® */(/)= V 0 eo *[("0- "oi-"o 2 ) ' + +9,] 

+ »,</>} 


(7) y (' ) = K 4 |a sin J(o»o - »*i - ( I + Af|) a» A2 ) / + 6 - 9 

+ (i + «,) e, +o 2 ] + />(/)| 

® z (/) = A-(p) y (/) 

®MO~ «/2 c..[^ + - VC °/ ( ^ ,</>] 

(lo) * 3 (/)~ 'Iz sin(A/|b) A2 t + 0 2 ) 

® *AGC (/) = ~ e G + ^{* cos [("0" < ‘ , AI _ ( l + *l)‘*'A2)' + °- 
+ (I +*,)#, + 0 2 ] +n (/)}. 

(j?) cU) - K c r{p) * AGC (/) 


Fig. 8-1 . The double-heterodyne receiver, with equations for signals at each point 


zone is noise; it has total power 2 (21/ n) 2 (1 — a 2 ). The 
limiter output signal- to -noise ratio pi is thus a function 
of a only: 

p« = m 

We shall investigate the way a behaves as a function of 
its input signal and noise a little later; for the present, 
let us express our answers in terms of a. 

The limitei output feeds a phase detector whose refer- 
ence input is rationally related to the free-running oscil- 
lator used to produce the second IF frequency; i.e., for 
some (relatively) fixed 0 2 , it may be taken as being of 
the form 

sin (Mi(i>» 2 t — Mi 0i — 0 2 ). 

(and Mi, like /, and f 2 , may take on a negative sign, if 
desired). 


Correspondingly, the detector output is 

1 /(f) — Kd^a sin { [«>o — w*i — (1 + Mi) <di, 2 ] t + 6 — 0 

+ (1 + Mi) 0i + 0 2 ) 4- n(f)^ . (8-2) 

The constant K d includes the gain of the detector, the 
limits ± l, etc. When the loop is in lock, the multiplied- 
up VCO output must have as its frequency 

<■>*! = <0(> — (1 + Ml) Uh2 • (8-3) 

The values of the IF frequencies are thus directly related 
tO fh2- 

fi = (l + Mi)f k2 

(8-4) 

f 2 = Mif h2 

which may vary with any fluctuation f k2 may have. 
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The heterodyne signal coming into the first mixer from 
the VCO through a frequency multiplier will then be 
one of the form 


{ 


x:(t) = cos \<ai n t + aMK (i Kveo 


EM 

p 


£sin (0 — 

MKycon v (t) \ 

P ) 


0 + (1 + Mi) 9\ + 0 2 ) + 


f] 


(8-5) 


including the equivalent VCO input noise, riv(t). The 
receiver estimate of the incoming phase function 0(f) thus 
satisfies the following relation: 


$ = «K rf Kv coM |sin[<9 - 0 + (1 + M,) 0, (t) + 0 2 ] 

+ ^J 

(8-6) 


If properly designed, the loop naturally adjusts itself, try- 
ing to zero the error; a linear analysis applies whenever 
the detector output error <£ = 0 — $+(1 + M,) 0, + 0 2 
is small; in this case, the loop estimate is 

a(t\ = aKjKycoMF(p) 

p + aK,K VC0 MF(p) 

x [*)+*?+ a + MM( t )+fc+^j^ y ]. 

(8-7) 


So when the receiver is tracking, = 0, — 0 4- (1 + M,)0, 
+ 0 2 is small, and the AGC detector works precisely as 
described in Chapter 7. Detection of phase modulation 
by a second synchronous detector also follows the same 
behavior exactly. Doppler-frequency extraction can also 
be implemented in a similar way. 

There are, however, some other minor differences be- 
tween the simple phase-locked loop and the double- 
heterodyne receiver, principally in the phase-detector 
outputs. It is necessary to distinguish between — 9 — 9, 
the loop tracking-error, and the quantity <f> — 0 — d 
+ (1 + M,)0, + 0 2 , the loop detector-error. The latter of 
these is the only one that need concern us if care is taken, 
as above, in proper information demodulation. In the 
simple no-heterodyne loops of Chapters 4, 5, and 6, <f> 
and were the same, but they are not the same quantity 
in the receiver above. 13 

Note that fluctuations in <j> occur (compare (8-7) with 
(5-24), for example) as if (1 + M,) 0, ( t ) were added to 
the input signal 0(f). Of course, the input signal itself has 
been generated by an oscillator whose output frequency 
has been multiplied to carrier frequency by some factor 
M-, and generally M-. » 1 + M,. Assuming that the 
oscillators associated with 0(f) and 0i(f) are of equal 
quality, insofar as phase purity is concerned, we observe 
that unwanted fluctuations due to (1 + M,) 0i(f) are 
masked by those of 0(f), the latter being roughly 
MJ{ 1 -1 M,) times as large. Hence, the contribution of 
0,(f) to <f> can usually be neglected. 


This response is the same form as that given in (5-2), 
except for the term (1 + + 0 2 , which contributes 

to the loop noise, and which may also cause a varying 
phase lag between 0 and 0, even in the absence of input 
noise. This term represents a heterodyne noise present in 
the loop. 

The fact that § does not follow 0 exactly (due to 
(1 + M,)0, + 0 2 ) need cause little concern so long as 
synchronous detection of information processes, such as 
doppler shift, phase modulation, etc., are also accom- 
plished using the same reference frequencies m 2 and 
M,<d *2 combined in the same manner as above. For ex- 
ample, the AGC detector uses a 90-deg shifted version 
of M ,<»a 2 to derive the stabilized amplitude of the signal; 
hence the input to the AGC filter and amplifier is 

fl<or(Q *1 

A J* 


The loop transfer function L(s) is exactly the same 
form as (5-3) with the parameter AK replaced by 

AK = aKdKvcoMF . (8-8) 

(Here we have made explicit inclusion of F, the loop- 
filter dc gain, previously assumed to be unity, in the 
event that some loss is present. In this way, we may 
assume that F(s) takes the previously treated form.) 


“There are ways of reducing the difference between <f> and For 
example, can be derived as «„td 4 - u h wherein a spectrally 
very pure standard reference frequency, such as one derived from 
an atomic-resonant oscillator, is mixed with a local free-running 
oscillator w 3 ; final phase detection then uses w 3 . The effect of this, 
in our terminology above, is that Afi is set to — 1, so that 0 
differs from e ^ only by $2 plus the standard oscillator phase, 
However, 0 a is a constant, and the standard is very stable; hence 
for most practical purposes, B and 0 differ by a constant phase lag, 
tunable if desired. 


— e G + AK d cos <f> 4- 
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8-B. Effects of Band-Pass Limiting 

In Chcpter 7 we discovered that AGC action may 
jitter the value of A by several db even at a given con- 
stant-input amplitude A(t). Without a limiter in the loop 
(A then replaces a in (8-8) above), the loop bandwidth 
would jitter in the same manner as A. But with a limiter, 
a is related directly to the predetection signal-to-noise 
ratio, generally a much more stable quantity. Further, a 
is usually much larger than A, so the loop gain need not 
be so great to give a desired bandwidth. 

Since the AGC alone is unable to maintain a constant 
receiver output level, limiting also tends to increase the 
effective dynamic range of the phase detector, and hence, 
that of the entire receiver. These reasons make it almost 
necessary to have IF limiting in a quality receiver. 



Fig. 8-2. Davenport's band-pass limiter zonal SNR curve 


1 


I 



f 


The limiter output zonal signal-to-noise ratio varies as 
a function of its input SNR in a known manner, derived 
by Davenport in 1953. However, Davenport’s result is a 
relation between the output signal strength and the total 
noise in the output zone. Generally speaking, a limiter 
tends to have a wider output noise bandwidth than does 
its input, whereas L(s) is usually chosen to have a much 
narrower response than either of these bandwidths. Neg- 
lecting internal oscillator noises and assuming that S n „(/<u) 
is fairly constant over the response region of the loop 
(i.e., w H » tv L ), one may derive a result similar to (5-8): 


^ _ S nn (0) Wl 

a- 


(8-9) 


The limiter output spectrum has some fiducial noise 
bandwidth Wf. 


approximation to this, which fits extremely well (see Fig. 
8-2) over the entire range of A 2 /N a tv H , is 


pi/pH 


0.7854 + 0.4768p« 
1 + 0.2384p„ 


(8-12) 


The corresponding signal amplitude suppression factor is 


a = 


0.7854p H + 0. 4768 /j £ *]'* 
1 + 1.024p w + 0.4768p ~ H 


(8-13) 


As the loop bandwidth w L is a function of loop gain, 
it is thereby a function of a. The bandwidth of a passive- 
integrator loop whose filter is that of (5-13) is related to 
the parameter r of (5-14). In the double-heterodyne re- 
ceiver with IF limiting, this r is 

f = (8-14) 

Tl 



1 

S„„ (0) 


(1 - a') 


(8-10) 


and therefore the loop fiducial bandwidth is approximately 


w L = 


aKjKv ccMF r* + ti 
2r 2 Ti 


(8-15) 


since, of course, the integral represents the total normal- 
ized zonal noise power in accordance with the implicit 
definition of n(t) in (8-2). Now as a result, 



Davenports result, which we mentioned earlier, is that 
pi ** 2 p H when the predetection SNR, p H = A7(N 0 w H ), 
is large, and pi*»(ir/4) p H when p w is very small. The 
actual formula for pi is very complicated. However, one 


At very high input SNR’s, a approaches unity, a result 
that causes the loop bandwidth to approach asymptoti- 
cally the value 


U)i( max) 


fCdKvcoMF T* + Tl 

2t 2 


j KdKvcoMF r 2 


foWLo 

ct 0 (r« + 1) ' 


(8-16) 


The phase variance can be related to input signal and 
noise parameters by introducing p H = A*/N 0 Wh into 


I 
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(8-11). With this done (again assuming negligible oscil- 
lator noise), the following expression for a 2 results: 


N n w L / Wh p H \ 

A 2 \wipi ) 

N 0 W L p 

A 2 


(8-17) 


The factor T = u> w p w /u> t p t is the limiter performance 
factor . One might at first presume, on the basis of Daven- 
port's result, that <r 2 would be degraded by a factor of 
it / 4 (about 1 db) at very low predetection SNR’s. But 
limiting tends to spread bandwidth, making w t > tv H . 
This effect tends to compensate somewhat for the w/4 
expected loss, and measurements indicate that, in truth, 
this is the case. 


As the predetection SNR becomes large, a different 
behavior results. There is a factor-of-2 improvement in 
<r 2 due merely to the Davenport phenomenon. But, in 
this case, it is not difficult to show that there is no 
further improvement, because the limiter’s output noise 
spectrum has the same shape as that at its point. This 
makes tvi asymptotically equal to w H . 


Going back to Davenport’s original analysis, Springett 
has computed the factor that we would have denoted 
here by (u>ipi/u/ H pH)- His result 14 , shown in Fig. 8-3, is 
well approximated by 


tv i pi _ 0.862 4~ 0.690p w 
W» pH 1 + 0.345p„ 


(8-18) 



Fig. 8-3. The rati * of band-pass limiter output signal- 
to-noiso spectral density to that at input. The 
reciprocal of this curve gives the limiter 
performance factor I\ 


The expression in (8-18) is the reciprocal of the limiter 
performance factor, T. 


The final resulting value of a 2 for a loop containing an 
IF limiter is therefore given by 


„ _ N 0 w l 


1.0 + 0.345 


\NqWh/ 


0.862 + 0.690 


G&)J 




NqWl 
A 2 ’ 


(8-19) 


“Springett’s result is based on the assumption that the predetection 
noise has an ideal band-pass characteristic. The results for an 
arbitrary unimodal noise density, however, to engineering accu- 
racy, are the same. 


This latter approximation is valid within % db over the 
usual operating range of signal levels ( p H < 1). More 
accurate analyses should, of course, proceed under the 
more exact expression. 
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CHAPTER 9 

BEHAVIOR OF PHASE-LOCKED LOOPS (NONLINEAR ANALYSIS) 


Most of the important characteristics of phase-locked 
devices can be predicted by the linear analysis presented 
in earlier chapters. But since the linear representation of 
a loop is predicated on high SNR’s in the loop passband, 
one can expect that, as this assumption begins to fail, 
significant accuracy is lost. At this point, it becomes 
apparent that another method of calculating loop per- 
formance is needed. 

There have been many who have tried to perform a 
rigorous, exact analysis of the noisy loop, but most of 
these have failed. Develet modified the linear model to 
a “quasi-linear” model, and Van Trees was able to approx- 
imate the loop character by a Volterra-series expansion. 
Both of these techniques are capable of extending loop 
analysis past the range of linear theory, but they ulti- 
mately fail when the loop gets sufficiently noisy. 

Tikhonov and Viterbi independently were able to solve 
for the exact steady-state behavior of the first-order loop 
by Fokker-Planck techniques. The second-order loop has 
been treated with some measure of success by Lindsey. 

This Chapter presents still another method, due to the 
author and applicable only to the calculation of output 
variances and spectra, subject only to rather general 
assumptions. One assumption necessary to the analysis is 
that the joint-probability function of the phase error can 
be expressed as a diagonalized orthogonal expansion. Such 
an assumption is known to be valid for Gaussian proc- 
esses, sine-wave processes, and the output process of a 
square-law device whose input is a narrow-band Gaus- 
sian process, and, in fact, it does not seem to be too 
restrictive. The second assumption is that the steady-state 
phase-error spectrum can be simply related to the spec- 
trum of the detected phase error. 


curacy improves as the loop-degree increases and as the 
loop filter becomes narrower-band. 


9- A. The Spectral Equation 

Since we are interested in the behavior of noisy loops, 
the small contribution of oscillator noise can be omitted 
from consideration. We focus our attention on the steady- 
state stationary equivalent phase process $, discussed in 
Section 4-B, and the process <j>. A slight rearrangement 
of the exact loop equation (4-4) gives 

<f> + AK sin 4 = 0 - K n(t) . (9-1) 

When we replace 4> by its stationary equivalent $, mul- 
tiply (9-1) evaluated at t = ti by its evaluation at t — t 2 , 
and average the resulting product, there results an equa- 
tion involving various correlation functions in 4>, sin 4>, 0, 
n(t), etc., whose double-sided Laplace transform gives 
the analytic continuation of the power spectral density: 

F($)F(-s) 

+ AK S*. *, n *($) F(~s) Bin ~~ s 

+ A 2 K 2 S 8 |n «in *(s) = ® 9e W + K fc S nn (s) . 


A result of Bussgang extended by Barrett and Lampard 
to a wide class of stochastic processes reveals 

R*. .m *(t) = i; R**(r) (9-3) 

where ?/ is a functional of p(4>)> given by 
1 / ,+00 

r) = — / (4> — E&) sin p(«fc) d<P . (9-4) 

® J -» 

We have used a 2 as the variance of 4». 


The analysis is not limited to Gaussian noise inputs or 
to wide-band processes, although such assumptions gen- 
erally make calculations easier. The Gaussian assumption, 
in fact, yields results that agree excellently well with 
the known exact theoretical behavior of the first-order 
loop over the entire useful range of signal and noise 
values. 

The advantages of this method are that it is simple in 
concept, and easy to perform and, further, that its ac- 


The condition on <t> needed to justify (9-3) is that its 
two-dimensional frequency function be expressible in the 
form 

p[*(f»). *(*»)! = PM* »)] plHtJ] i*4. [*(*i)] /•[♦(**)] 

n-0 

(9-5) 

in terms of polynomials /„($) orthonormal with respect to 
the kernel p(4>). When such an expansion is valid, 

S*. .m •(«) = S*. ,m ♦(-«) = ijS**(s). (9-6) 
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Even if (9-5) is not valid, the only modification needed in the theory is that i> in 
(9-6) must be replaced by 7 (s). 

Next, when <t> is passed through the nonlinearity, sin <t>, there is a distortion 
of the zonal spectrum from S*« into S, ln * *(s). Hence we define a function 

y (s) y(— s), which relates the way this occurs: 

ym - t) = . ( 9 - 7 , 

Substituting these latter relations into the spectral equation, we obtain an expres- 
sion for the spectrum of namely 

= ’""S 2 Se* (s) + K 2 F ( S ) F( — s) S Mn (s) 

K ' -s ! + r,AK [sF (-s)-sF (s)] + A-K 2 F(s) F(-s) y (s) y( -s) (9-8) 

from which the variance of $ can be computed: 



Generally (9-9) will be transcendental, since both 77 and y(s) are nontrivial func- 
tions of a 2 . 


If y(s) y( — s ) is reasonably well-behaved, and if F(s) is a very narrow-band 
device, then | y(/«i>)F(/w) | 2 y 2 (0) | F(j<o) | 2 . Similarly, under the same condition, 

S n n(/oj) | F(;'u)) j 2 «= N 0 1 F(/'io) | 2 , where N 0 = S„„( 0). These relations, then, are 
particularly applicable when F(s) is a filter with a high cutoff rate. In fact, as 
long as the noise n(f) is wider-band IS than F(s), the system performance is essen- 
tially the same as if the noise were white, and the approximate S*« in this case 
is given by 


S**(s) = 


~s 2 S„ (s) + K 2 N 0 F(s)F(- s) 

-«* + vAK [sF(-s) - sF (»)] + A 2 K-’y 2 F(s) F(-s) ' 


(9-10) 


This method of 


analysis is thus called a linear spectral approximation. 


As usual, there are two terms that contribute to the phase error: one due to 
the input signal component, and one due to noise. The term we wish to explore 
more fully here is that due to noise: 


/ 


C (s)= K'-N n F(s) F(—s) 

’ -s- + V AK [sF ( -s) - sF (*)] + (AKy) 2 F(s) F{-s) ' (9 * U) 

Concerning the spectral approximations, one further point can be made: In the 
linear and quasi-linear models of the loop, <ft is approximated by a linear function 
of 4>, so in turn, S iln # *(s) and S* „ „*(s) are proportional to S**(s). Hence, the 

method above yields the linear results when the value for 17 and y(j) are replaced 
by the common value given them by either the linear or the quasi-linear model. 
However, while including the linear cases, the linear spectral approximation is 
a more general concept, because 17 and y are not necessarily equal. It is thus 
more apt to provide results closer to the exact behavior of the loop. 

"It must be remembered, however, that n(t) is the result of multiplying n, (t) by cos (**+•), 
and hence will be somewhat wider-band than n ( (t). 
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9-B. Calculation of r/ 

The calculation of q as the constant relating S*« to 
S* in * involves first knowing p(4>), and second, knowing 
that p[$(ti), $(# 2 )] has the diagonal form of (9-6). 


For the density in (9-12), the value of y can be readily 
computed: 

v = e““ 2/2 . (9-14) 


At very low noise levels, assuming a Gaussian channel, 
the density for p(<£) certainly approaches a Gaussian 
density, as does p(4>), and all analysis methods give com- 
patible results; at high noise levels the density on <\> (mod 
2 tt), viz., p(<f>), approaches a uniform density of l/2w on 
the interval ( — «•). In what follows, we denote <f> 

(mod 2ir) merely by <f>. 

As an approximation to this behavior, for definiteness, 
we may assume that $ is Gaussian, with variance a 2 ; such 
an assumption satisfies both limiting conditions on the 
density of 4>, which is then given by 

I I +00 

JtSr.?.“ p[_( ^ + 2m ' )y2j ’ h 

for | | < ir 

0; for 1 4> | > ir . (9-12) 

The variance <r 2 of the ^-process is related to the variance 
a 2 by straightforward evaluation. The lesult is graphically 
presented in Fig. 9-1 along with the approximate formula 

<r = y j L ~ exp (1 + 0.13o 2 )jj . (9-13) 


9*C. Calculation of y 2 

For loops in which F(s) is very narrow, y(s)y(— s) can 
be replaced by its value at s = 0, viz., y 2 = y 2 (0). In the 
case that F(s) has a very wide bandwidth, we can still 
use (9-11), but a better approximation to y(s)y(— s), in 
this case, is the ratio of the total power in sin 4> to that in 
4>, rather than y 2 , the ratio of powers around zero fre- 
quency. This total power value, y 2 p , is 


5 _ R«ln • »in<(0) 
y ‘> 


(9-15) 


Since we have assumed that $ is a Gaussian process, 
Price’s theorem can be used to compute the autocorrela- 
tion function of sin $: 


R.m *(r) = e - ®* sinh R*(r) . (9-16) 


The relation above is certainly valid when o 2 is small, 
because we are sure, in this case, that 4> is nearly Gaus- 
sian. For larger a 2 , (9-16) is only approximate. But insofar 
as the Gaussian assumption is justified. 


y* 


e-“ 2 



sinh R«(r) dr 


fl*(r) dr 


(9-17) 



Fig. 9-1. Relationship between the variance of ^ (mod 2») and that of a Gaussian stationary process * 
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for a loop with a narrow-band F(s), and 

y% = a ' 2 e_ “ 2 sinh a 2 (9*18) 

for a loop with a wide-band F(s). 

The former quantity, y 2 , cannot be evaluated directly, 
since R*(t) is not known until y 2 is, and vice versa. How- 
ever, if loop noise is not completely overpowering, y 2 is 
probably not as sensitive to the shape of R»(r) as it 
is to the value of a 2 . As a verification of this conjecture, 
the several different forms of R*(r) shown in Fig. 9-2 
produce values of y that are in very close agreement 
over the range of interest. The upper curve is equivalent 
to y tp . The values of y are probably very well approxi- 
mated by that in the middle: 


1. Linear Loop 

Under the usual Unear assumption, y tp = 1, so ti>£ ((q) 
reduces to w L , the linear loop bandwidth, and, of course, 
the usual answer (5-8) results. 

2. Quasi-linear Loop 

When the constant gain A in the linear loop is replaced 
by the equivalent statistical gain Ae~ a,/2 derived by 
Develet, i.e., when y J = e - ®’, the result is 

°* = IF (9*23) 

Notice that there is no solution for a* here when 
N 0 K/2A > 1/2. For N 0 K/2A < 1/2, there is an equiva- 
lent bandwidth 

AK 

= X e- a/2 . (9-24) 

3. Spectral Approximation 

If we use the analysis presented above with y tp as 
given in (9-18), the phase variance of the equivalent sta- 
tionary process 9 is related to the Unear-loop prediction by 

aer** 1 * (sinh a 2 )* = ^ . (9-25) 

Here there is a solution a* for every value of N«X/2A. 

4. Reduction Modulo %w 

The real indication of system performance, as we have 
agreed, is given not by the parameter o* (the variance of 





Fig. 9-2. Variation of tho parameter y as a function of 
•.he Gausstan variance a 2 , for various forms of Ih (t). 
Note that, fora 2 < 1, there is not a significant 
dependence on tho form of JU* It). 


y = 2 <r 2 e-® 2 ' 2 sinh I 


(1 - e*“ s ) . (9-19) 


To see how well the spectral approximation method 
works, let us calibrate it by a few examples. 


9-D. Behavior of the First -Order Loop 

The spectrum of $ for a fir^^order loop is approxi- 
mately \ 

\7 K' 2 

O /ft 


S** = 


— S 2 -f (AKy tp j' 


(9-20) 


for we merely insert F(s) = 1 into (9-11). T Ve use y t P to 
approximate y(s), because F(s) is wide-band.\Integration 
of (9-24) yields the estimate v 


No K Np / Ay„K \ 

2Ay <p A 2 y 2 p V 2 )• 


This equation compares in form with (5-8), except tha\ 
the bandwidth of the loop has apparently changed to a ' 
new equivalent value because A reduced to Ay tp : 




_ Ay >g K 


(9-22) 
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<I>), but by a 3 , the variance of 4> (mod 2m). These are 
related as shown in Fig. 9-1. Each of the various per- 
formance approximations, reduced mod 2m, is compared 
with the actual loop behavior (Fokker-Planck method) 
in Fig. 9-3; notice the closeness with which the spectral 
method predicts the actual behavior. We have also indi- 
cated that the spectral approximation method should 
yield results that improve as the filter F(s) becomes 
higher-order and narrower in bandwidth. Thus, the first- 
order system should be the one predicted with least 
fidelity — but, as shown in Fig. 9-3, even this proves to 
be very good. 

The conclusion to be reached at this point, then, is 
that the linear spectral approximation method provides 
a very accurate means of predicting loop phase devia- 
tions (mod 2m) over the entire useful range of input signal 
and noise values. In addition, we have shown that the 
linear analysis is valid whenever A i /N 0 w L > 10. 

\ 

X V-E. Calculation of tha Behavior of the 
\ Second-Order Loop 

\ 

VXe now show the performance of a phase-locked loop 
whost^loop filter is the passive integrator of (5-13). This 



Fig. 9-3. Comparison of linear, qua si-lli' oar, and linear- 
spectral approximate methods with the actual 
behavior of the first-order lot p 


filter function inserted into the spectral equation (9-11) 
produces the density 

S**(s) - 

K*N„( l-r*«*) 

7*? - s*"tl"+"2AK,(r, - ~) + (AKy r 2 ) 2 ] + (AKy) 1 

(9-26) 

Ihe density (9-26) is quite different from what one ob- 
tains from a simple linear model of the loop. Still, we 
caB 'iise the same types of ideas to express the output 
phase variance. For one thing, S«*(/n>) has a fiducial 
bandwidth 


( eq ) 


1 

2rS**(0) 


/ 


■too 

S44 ( /<!)) dot . 

40 


(9-27) 


Note that this is computed in the same way as w L , except 
that we must use S**(/«>) as given in (9-26) rather than 
| L(/») j 3 , as one would do if the linear model were in 
effect. In these terms 


_ N„w L(tq) 

~ My)> ' 


(9-28) 


Again, it is evident that there are two effects in the 
second-order loop that deviate a * from the value it would 
take in a linear analysis: w LlHI ) must be used rather than 
Wi, and A is reduced to Ay. 


The value of u> can be found by integrating (9-27): 

(1 + ry) 

n /i i- > rj \ T i r 4y - q)(i - r*A>)"P 

2r, (1 + r,/fr iy ) 1 + v<1 + ^ y J, _ 




1 + fy 

L r r J 


(9-29) 


The approximation follows the natural assumption 

ft iy » r,. 


A “damping factor 1 * can be defined almost in the same 
way as it is for the linear loops, a way that reduces to 
the linear care for high input signal-to-noise ratios. We 
have already derived S««; it is of the form 

«"W - 7^T-c ■ 


The roots of the denominator govern the “damping be- 
havior” of *(t). If these roots are r,, r„ -r tl — r,, then 
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we can define, for the left-half-plane factor s 2 -f (r, + r 2 ) s 
+ r,r 2) an equivalent damping coefficient 

We relate r,r 2 and (r, + r 2 ) to the parameters of S #*(«) by 
noting 

(9-32) 

(fi + r 2 ) 2 — 2f!r 2 = T" 2 [1 + 2AKt)(t2 — t,) + (AKyr 2 ) 2 ] . 


3. Linear Spectral Analysis 

In a still more general analysis, we set q = e~“ ,/2 and 
Y = a' 2 (1 — e - ® 2 ). The loop behavior is then approxi- 
mated very closely by 

NqWl _ (1 — g ~°*) 2 W L 

A 2 a 2 w L(tQ) 

w L(tq) 1 + r(l - e^)/a 2 

~(l + ,)[l + 2 51ii^^2]“ 

(9-35) 

U (l-«TTi , 2a 2 (1 — e~* — a 2 e~ a>/2 )~\ * 

K ~ a L J * 


Hence, the equivalent output damping factor of the loop is 

[d + Vr„,)- + ^ ^ ~ ’MJ 

+ (M3 , 


Again, the same type of behavior results: and 

are less than their linearly predicted value, while a 2 is 
greater. Although the quantities in (9-35) are generally 
transcendental in a 2 , nevertheless Newton’s rule is easily 
applied by a digital computer for solutions. Figure 9-4 
relates various parameters above to corresponding values 
of a 2 , based on a constant value of r (viz., r = 2). Figure 
9-5 then shows how a 2 compares with the value pre- 


1. Linear Analysis 

It is easy to see when q = y ~ 1 (the linear analysi* 
parameters) that (9-28), (9-29), and (9-33) reduce to (5-8), 
(5-18), and (5-20), respectively. 


2. Quasi-linear Analysis 

The quasi-linear analysis replaces A everywhere by 
Ae - ® 2 2 (i.e., y= q — e*® 2/2 ). Thus the behavior predicted 
in this case is given as 


_ W 0 tOi,<eq) 

ae -g— 


tt>£(eq) 


1 + r e-®v 2 

2t 2 


(9-34) 





From this it is important to note that while Landwidths 
and damping are less than the value predicted by die 
linear model, the overall value of a 2 is larger than its 
linearly predicted value. 



Fig. 9-4. Comparison of linear and nonlinear theories 
for second-order, constant linear bandwidth loop, i.e., 
the value of r is kept constant at r 0 = 2. 
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dieted by a linear analysis as a function of signal power 
for a fixed N 0 ; since r varies with input signal level, we 
have arbitrarily 16 chosen a value r 0 = 2 at AJ =3N 0 /2r«. 


"Reasons for this choice are given in Chapter 10. 



Fig. 9-5. Comparison of phase-noise variances by linear 
and linear-spectral approximations. The noise density 
is fixed, and the signal level is varied. The value of r 
is taken as r 0 = 2 at a reference signal level of A 2 
= 3N„/2t 2 = N„w to . Note that the ultimate roll-off 
is 5db/decade, rather thou 10 db/decade, as in 
Fig. 9-4. Note also that even the linear approxi- 
mation produces some curvature of <r 2 near A 2 . 


Figure 9-6 illustrates the way w t ., ( and t, 

vary; parameters here are the same as those in Fig. 9-5. 

4. Conclusions About the Second-Order Loop 

The behavior depicted in Figs. 9-4, 9-5, and 9-6 indi- 
cates that a linear theory can be expected to yield satis- 
factory accuracy whenever A 2 /N 0 w L > 10. Beyond this 
point, the linear spectral approximation probably agrees 
with actual results, if we may extrapolate the results 
obtained for the first-order V op, for o 2 < 1.5 rad 2 . This 
figure lies beyond the useful ange of most receivers. 



Fig. 9-6. Variation in bandwidth and damping param- 
eters as a function of signal strength. The value of r at 
a reference signal level A 2 = 3N„/2r 2 = N 0 w,, 0 
was taken as r n = 2 . 
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CHAPTER 10 

DESIGNING A DOUBLE-HETERODYNE TRACKING LOOP 


So far in this work, we have concerned ourselves, for 
the most part, in predicting the performance of a receiver 
when its parameters were specified. Exceptions to this 
dealt only with ways of choosing some of the receiver 
parameters, given the others. We may now pick a set 
of parameters, and, through the formulas that have been 
presented, predict quite accurately how the loop will 
operate. Once the receiver is built, it performs just so. 

What we require for design, however, is an effective 
method for picking a nominal set of values that produces 
a desirable tracking function over a reasonable variation 
of the input signal power. 

As long as the operation of the loop is to remain within 
the linear-analysis region, the task is much simplified, 
because there have been analyses put forth for extracting 
an optimum result from the linear loop. Nonlinear op- 
timization in most cases is difficult, if not analytically 
impossible. 

It is customary to optimize, insofar as possible, receiver 
performance at the “worst-case” parameter values; in this 
way one is sure that, while perhaps not operating optimally 
at any other set of values, the loop will do the best that 
it possibly can in those cases that require it most. This 
procedure is purely an arbitrary one, in that it is apt to 
change in accordance with the philosophy of the design 
engineer, and with the particular mission for which the 
receiver is intended. 

What design rules are given here will be somewhat 
abstract, but perfectly general, so that the designer may 
issue his own philosophy in choosing values. 

JO- A. Definition of Receiver Threshold 

The words "receiver threshold” conjure up a different 
image to each engineer: to one, it is that shaky, ill- 
defined signal level at which the receiver transits from 
operability to nonoperability; to another, it is some point 
at which the receiver operates with the least acceptable 
reliability; to still another, it may mean the signal level 
at which the receiver exits from its linear range. There 
are no absolutes when it comes to defining such a point; 
in r ;h of the cases above, the engineer meant that the 
system hit “threshold” when its performance was in some 


way no longer tolerable, when the receiver no longer 
produced meaningful data. The trouble with such defini- 
tions tied to performance is that they tend to be sub- 
jective — what is tolerable in one case is not in another. 

The word “threshold,” as we shall use it here, is a 
precisely defined quantity, which can be subjectively 
interpreted as desired. Some give it the name “design- 
point threshold” and, in some cases, it has been called the 
“absolute receiver threshold.” The concatenation of modi- 
fiers does not seem to be necessary at all, and we shall 
use the word “threshold” in only one sense. 

Specifically, the loop threshold is defined as that input 
signal power A* given by 

_ N 0 (r 0 + 1) 



= N 0 w Lo . ( 10 - 1 ) 

Here, w Lf> is the value of linear-loop fiducial bandwidth 
at threshold, i.e., computed by (5-18) using the threshold 
value of r, r 0 . One must not be mistakenly led to believe 
that the loop is operating linearly; is merely the 
bandwidth a loop would have if A 0 were to lie in a linear- 
analysis region. In actuality, A 0 is a signal level at which 
the linear theory does not apply, for, in fact, <r* > 1. 

In the notation that is to follow, all quantities sub- 
scripted by “0” refer to the value of that particular 
parameter at threshold. 

70-B. Tracking Loop Performance of the 
Double-Heterodyne Receiver 

Once the values of the receiver’s physical parameters 
KdKvcoMF, n, t 2 , and W H have been determined, per- 
formance proceeds as specified by the formulas 17 in 
Chapter 8. But with a fixed receiver, the loop bandwidth 
changes as a function of signal or noise level. Rather than 

"Although Wh is a measured physical quantity, it may not be 
equal to Wu if the loop is not operating at the frequency pro- 
ducing H(iu) Care must often be taken to assure that H(s) 
is properly tuned for maxi>num response at the loop operating 
frequency. 
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calculate each receiver characteristic separately, it is 
easier to evaluate the performance in terms of a reference 
signal level and to relate the remainder of the loop be- 
havior to this point. One convenient reference point is the 
threshold signal level. Specifically, we define the ratio 


m = 



( 10 - 2 ) 


as the loop performance margin. In terms of this m. 


Ph = 


A 2 

N„tc„ 



(10-3) 


As a result, all of the 1«* p characteristics are expressible 
in terms of m, w La , w„, and tJt x (and this latter quantity 
is not needed if t, » r,). Thus a design consists really of 
specifying these three (or four) quantities. 


However, it is often necessary to compute the perform- 
ance of a receiver when tc H , KdKvcoMF, t 2 , and r x are 
given — quantities that can be measured. Most of the 
formulas we have given are in terms of r s, w t ‘s, etc. How 
does one proceed with the set of parameters above? We 
merely need to find r 0 , and the rest then follow. When 
10 u>t 0 < w H , a 2 is very well approximated by ww Lo /4w H ; 
from this, a slight manipulation of (8-14) produces the 
desired result: 


17 = 


*(K d K VV oMFYr\ 


8 T* U)„ 




We can also express U in terms of r 0 , 


(•*•0-4) 


TJ _ fo(r 0 + r 2 / Ti) 

r„ + 1 


to 4- 1 


(10-5) 


70-C. Nonlinear Behavior of the 
Double-Heterodyne Receiver 

The formulas in preceding chapters can be made to 
apply to loops with predetection bandpass limiters by 
making a few minor alterations, and we have seen how 
these changes come about in Chapter 8: first, AK is re- 
placed by aK d KvcoMF, as in (8-8); and second, the ex- 
pression (9-28) must be multiplied by the limiter 
performance factor T = w H p H /wipi, shown in Fig. 8-3. 


These steps lead to a set of equations which characterize 


the double-heterodyne receiver: 


<‘- = m fe) =m '"" 

(10-6a) 

[ 0.7854p„ + 0.4768p 2 -|* 

a l 1 + 1.024p„ + 0.4768p 2 J 

(10-6b) / 

1 + 0.345p„ 
0.862 + 0.69 0p„ 

(10-6c) 

aKdKvcoMF(0)rl a 

r = = — r 0 = t‘B 2 

ti a 0 2f ^ 



(10-6d) 


= Y {l - ex p[^ (1 + 0.13a 2 ) ] j. 

~ a 2 (10-6e) 

y = (1 - e - ® 2 )/a 2 (10-6f) 

v = e-" 2 '' 2 (10-6g) 


w L = 


1 + r 1 + r 

2 '-( 1 + m) 2 ” 


1 + 




L 1 + r 0 J 


u>h 0 


(10-6h) 


U>Z.<e q) — 


2r 2 fl + — ^ 

~1 1 2 (? ” 3)0- “ T */r 0~ 

V4 

V »W 

y 2 r(l + — V 



L V J 



(10-6i) 




£ "-III 
2 ~ 2 


(10-6j) 


= y 1+ n + . % - - agan l* 

2 V ryn) y* L f y 2 ( 1 + T */ r m ) 2 Jj 


*=* £ (*>)“ — i v* 


(10-6k) 
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T 2 


(10-61) 


The subscript “1” on T, y, -q, etc., refers to evaluation at 
a 2 = 1. Under the usual assumption y,r„Ti * t 2 , the approxi- 
mate value for m, is 


2 N I i‘q 1 U _ 1 _ / WLteqj] Q\ 

AY - m\ w L( , y- ) 


(10-6m) 


2,w L _ / n V* 2ic, m 
r + 1 \u» / r„ f 1 


(10-6n) 


The use of y as given in (10-6f) requires » r 2 ; otherwise, 
y, p of (9-18) should be substituted in its place. Notice that 
the only quantities needed to specify everything in (10-6) 
are the margin m, the predetection SNR p H( , = r,„ 

and to/tj. 


m t 


-[whr]"{ 1 + [ 1 + 


4(ro+l) 1 *\* 

r.rj J / • 


( 10 - 8 ) 


The value of a"- is within 0.5% of unity when a 2 = 1; 
hence to the accuracy we need, the values of mi and r, 
above give a 2 — 1. The curves in Fig. 10-1 show how mi 
varies as a function of r„ and t 2 /t x . As nears unity, it is 
necessary to replace y by y tp (see 9-18). The values given 
in the Figure result b' - assuming that y, in (10-3) takes 
the value 


The performance observed here is very much the same 
as that exhibited in loops without limiters, except that the 
value of KrfKvcoMF required is drastically less than the 
value of K required when there is no limiter in the loop. 
This, of course, is due to the fact that while A n is very small 
(the threshold rms signal level) ct n is many orders of magni- 
tude greater, approximately equal to the predetection SNR. 

70-D. The Signal Level Producing cr 2 = 1 

Since it is doubtful that one would ever expect a loop 
to operate usefully with a value of a 2 greater than unity, 
the specific value of a\ is more a point of academic interest 
than anything else. Furthermore, if we may judge from 
the first-order loop result, the threshold value of phase 
noise predicted by the linear spectral theory would 
probably not be a very accurate one, for the theory 
begins to fail somewhere near this point. 

Holding the rest of the loop parameters constant, we 
can solve for the value of margin (call it m,) that produces 
a 2 = a\ =1. This is a more meaningful quantity, and 
certainly much more accurately calculable by the methods 
in Chapter 9. It follows from (10-6) that when a 2 = 1, the 
value of A,/A 0 = mj 4 satisfies the equation 


= 0.6321 ^1 - ^ + 0.60653 ^ . 

(10-9) 



Fig. 10-1. Variation in margin producing a 2 = 1 a* a 
function of threshold design parameter, 

to = a.oKdKvcoMFr \/ t\ 


A± _ 
A„ 


r,(r 0 + r 2 /Ti)(A 1 /A 0 + 1/foyi) 


yi(r 0 + l^A/Ao + T 2 /r 0 yiTi) 1 + 




2(yi ~~ iyi)(l Tg / ti)(Ai/ Aq) 

y\ r 0 (Ai/A 0 + Ti/ySoT^j 2 J 


Vi 




a, At 

f i = — fo -7- f o • 

Oto -Aq 


(10-7) 
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70-E. Choice of Receiver Parameters 

Choosing values of N,„ w Lo , r n , K^KvcoMF, and w H speci- 
fies a typical receiver design. Values for n, t 2 , p H(l , oto, and 
T„ follow directly from (10-1) and the first four parts of 
(10-0). The value of r 2 /r t is a measure of how well the 
imperfect integrating loop filter performs in the loop; with 
t 2 /ti = 1, for example, the loop filter is merely F(s) = 1, 
and all the results in (10-6) reduce to the first-order loop 
equations. On the other hand, the results in (10-6) are 
approximately the same for all values of less than 
about r„/10. There is some flexibility in choosing loop gain 
parameters because K d , K V co, M, and F appear as a single 
factor in the theory. 

However, the only things an engineer needs to know 
in order to know how well any receiver will perform are 
m, r 0 , p Ho , and perhaps t 2 /t,. This latter set of quantifier 
depends on several things: r 0 sets the loop damping factor, 
to be chosen in accordance with the discussion in Section 
6-C4; the value of p Wo depends on w Lo and w H ; w Lo is 
determined by considering the desired acquisition ease, 
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Fig. 10*2. Comparison of linoar and nonlinoar approxi- 
mations to loop rms phaso orror, as a function 
of loop margin 


the doppler tracking capability, the state of the art in VCO 
design, the spectral purity of the received carrier, the 
expected incoming A 2 /N 0 , and the amount of phase jitter 
that can be tolerated at that operating level. In addition, 
T\ should be chosen (Section 6-C3) so that the phase error 
due to the expected doppler rate is increasing slowly 
enough that only infrequent retunings of the VCO are 
required. 

Once r„ and p H() (and perhaps t 2 /ti) are given, the re- 
ceiver performance in terms of margin can be found from 
( 10-6); a typical plot of a vs m is given in Fig. 10-2, wherein 
linear approximations are compared with actual behavior; 
Fig. 10-3 shows how bandwidth and damping factor 
change as a function of margin. With assumed values of 
N 0 and w L(j , threshold is specified (10-1), so the curves are 
easily referenced to the actual incoming signal power. 



0 3 10 IS 20 25 SO 38 

MARGIN m, db 

Fig. 10-3. Variation of loop bandwidth* and damping 
factors, as a function of loop margin 
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PART II 

SUMMARY OF PHASE-LOCKED LOOP DESIGN CONCEPTS 

In this second part, the important definitions, concepts, and formulas are col- 
lected chapter by chapter. The notation used in these formulas is listed in the 
Appendix with names, units, and text references. Many of the formulas here are 
only approximations of more accurate ones in Part I. 

Since not all of the subdivisions in Part I are referred to in Part II, there are 
some discontinuities both in the headings and in the figure and equation numbers 
of Part II. 
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CHAPTER 1 

INTRODUCTION AND HISTORY OF THE PHASE-LOCKED LOOP 


The first serious application of the phase-lock principle was as a synchronizing 
device for television in the late 1940’s. Since then it has evolved as the heart of the 
most sensitive, versatile, and flexible receiver in existence today. There are many 
approximate analyses for predicting the behavior of the loop (Fig. II-l-l); some 
of these are set forth in the ensuing chapters. 



Fig. 11-1-1. Basic configuration of a simple phase-locked 
loop. The mixer output, filtered by FIs), is used to 
control the frequency of the voltage- 
controlled oscillator (VCO). 








FUNDAMENTAL CONCEPTS 


2 -A. Statistics 

Process: a time function, x(t). 

Deterministic process: a process x(t) for which x(t) is 
specified at every value of t. 

Stochastic (or random ) process: a set, or ensemble, of time 
functions (*(*)}, any particular member of which is 
observed according to some probability law. 

Stationary process: a stochastic process whose statistical 
behavior does not depend on its time origin. 


Cross-spectral density S^d <o): similarly defined as the 
Fourier transform of flr»(r). 

Mean value of x(t): 

mW = E [*(0] 

= n (if x(t) is stationary). (2-10) 

Variance of x(t)-. 

*l(t) = E[x(t) - M (t)]* 

= E[x*(t)]-S(t). (2-11) 


Ergodic process: a stochastic process wherein time aver- 
ages produce the same results as statistical averages. 

Time-autocorrelation of the function x(t): given a function 
x(t) 

31**(t) — lim f x(t)x(t+r) dt. (2-3) 

T -+00 Ai J_ T 

Statistical autocorrelation of a stochastic process (x(t)}: 
BJfiut,) = £[*(#.)*(<,)] (2-4) 

= R**(t, — t t ) (if x(t) is stationary) (2-5) 


2-B. Linear Filtering 

Unit impulse response of filter H (Fig. II-2-1): The output 
h(t) of H when the input is i(t), the unit-weight Dirac 
impulse function (see Fig. II -2-2). 

Linear filter: any filter whose input x(t) and output y(t) are 
related by 


y(t) --- J x(t x )h(t - tjdti . 


(2-14) 



= l TLAti - t 2 ) (if x(t) is ergodic). 


(2-6) 


Fig. 11-2-1. Filtering device 


Cross-correlation R n (ti—t 2 ): obtained by replacing x(t 3 ) 
by y(t-) above. 


Spectral density of a function x(t): 

3Ub) = j" HUder^dr. 
Spectral density of a stationary process {*(*)}: 


SJfm) = f" R»»(t) e-i^dr. 


(2-7) 


(2-6) 



Fig. 11-2-2. iHponu of filter to a unit-impulse function 
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When x(t), y(t), and h(t) have Laplace transforms X(s), 
Y(s), and H(s), respectively, linearity reduces to 

Y(s) = X(s) H(s) . (2-15) 

When x(t) is a stationary process, 

SM = \H(1»)\*S„(j»). (2-17) 

2-C. Noise Bandwidth 

White noise: a random process n(t) such that 

Rnn(r) = N„ 8(r) (volts*) (2-20) 

Sn»(j< ») = N 0 for all — oo <«< +oo (volts Vcps). 

(2-21) 

Effective noise bandwidth of a linear filter H(s): 

•fcf" |tf(/<,)|*da, 

W„ = 2B„ = (‘T»s) ■ (?-24) 

Note that W„ is the two-sided bandwidth find B H *s the 
single-sided bandwidth of H(s) in cps (me Figs. 11-2-3, 
II-2-4, and II-2-5). The output noise power with a white 
noise input is 




2-0. Sinusoidal Filter Inputs 

If 2 1 ' 4 P sin(<i> 0 1 + 0) + n(t) is put into H(s), where n(t) is 
white with spectra! density N 0 , the output SNR of H(s) is 


Ph = 


nW,, 


|H(/ W0 )|* _ P 

|»(«|L " ' 


(2-35) 


N = N,W„ |H(h)| L. (volts*) . (2-25) 



2-f. Fiducial Bandwidth 

Fiducial bandwidth: equivalent (2-sided) bandwidth re- 
ferred to a s;-. > if!c frequency «■><>, such as the carrier 
above, rather than that « producing |H(/«)| , 


w H — 2 bn 



(cps). 

(2-37) 


It is this fiducial bandwidth that is needed in most of 
phase-locked loop theory, rather than the equivalent 
noise bandwidth W H . 



2-F. Band-Pass Mixers 

Heterodyning (Fig. II-2-6) does not affect signal-to-noise 
ratios so long as the heterodyning frequency and the 
carrier frequency are independent. However, when the 
two are the same, the output signal-to-noise ratio p is 


P = 


P 

N'Wd 


(2-52) 
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Fig. II -2-4. Double-sided frequency response 
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•Ji. co»lw A t+8 3 ) 



where P = input 'Ignal power (watts), N 0 = 2-sided input 
noise density (watts/cps), w D is the fiducial bandwidth of 
the detector (cps), and ^ is the phase angle between the 
carrier and the detecting sinusoid. The term E( cos 2 <f>) is 
the coherence factor of the detector. Image noise must be 
eliminated by post-heterodyne filtering (Fig. II-2-7). 


2-H. Molt* 

Johnson (or thermal) noise: noise due to thermal agitation 
of electrons in a .resistor. There is a noise voltage t;,(t) in 
series (see Fig. II-2-8) with each such resistance R, with 



S VnVn (j») = 2kTR (voltsVcps) (2-62) 

where k is Boltzmann’s constant and T is the resistor 
temperature in degrees Kelvin. This is Nyquisfs Law. 
Maximum noise power obtainable from such a source 
in a bandwidth W v is 

N = W.v =■ I TBs (watts). (2-6ci) 

(Compare this with (2-25), for a filter with 1 -ndwidth 
Wjy and gain |H(/«>)| J,„.) 

Hence for a white thermal noise the matched-load noise 
voltage appearing across the load resistance R has uni- 
form spectral density .V 0 = 1 kTR voltsVcps. 







S tt (y«)« 2*r<? 


Fig. 11-2-7. Choosing filter bandwldths to avoid th« imago 
noiso problem 


Fig. ll-2-t. Thevonln and Norton equivalent circuits of 
noisy resistors 
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CHAPTER 3 

FORMULATION OF THE LOOP EQUATION AND 
BEHAVIOR IN THE ABSENCE OF NOISE 


3- A. The Basie Integro-Differential Equation 

The basic equation relating to the device in Fig. II-3-1 
is 

A(i) — + AK ^ sin (3*9) 

where A = rms signal amplitude, volts 

K — total open-loop gain, sec* 1 volt -1 
0(t) = input phase (unction, rad 
<f>(t) = phase-error function, rad 
p = d/dt, the Heaviside operator. 


3-B. Tracking When the Loop Error Is Small 

If |4>| < w/6, then sin <f>(t) =« so 

m~ p+ A ^ p - *«) . (3-io) 

V 

Considering 0(t) = d(t) radians, i.e., a doppler-shift func- 
tion only, the steady-state phase-error is 


<f>„ = lim 
8-* 0 


s*D(s) 

s+AKF(s) 


(3-12) 


where D($) = _C[d(t)]. If d(t) - and if F(s) has 

l poles at the origin, as F(s) = q($)/s l p(s), then for n > l + 1 

f n-i—l j i (W6) 

For n < l + 1, = 0, so perfect doppler tracking ulti- 

mately results. When n = l + 1, there is a constant steady- 
state error 




- AnP(O) 
AKq( 0) 


(rad). 


(3-17) 


3-C. Acquiring Lock in the First-Order Loop 

A phase-plane diagram of first-order (F(s) = 1) lock-in 
is given in Fig. II-3-2. Lock-in proceeds as 

(rad/sec) = AK = 2 W L (cps) (3-23) 


*..= = sin_1 (w:) M> ( 3 - 24 ) •’ 


t»ca ^ 


In 


W i COS <f>gs Sjoclf 


(sec) (2-37) 


where ft mit is the maximum value of (l 0 , the initial loop 
frequency offset in rad/sec :'or which the loop ultimately 
locks, facq is the maximum time in seconds required to 
achieve lock within S| 0C k radians, and W L is the 2-sided 
loop bandwidth in cps. 


3-D. Acquiring Lock in the Second-Order Loop 
with Passive Loop Filter 

With a loop filter of the form 


"W = TT^ ^ 

where r,, r 2 are in seconds, for d(t) = 9 0 + (lot, there 
results 

= «in" (^) = sin- (^-) (rad) . 

(3-30) 

But for lock-in to occur for all initial conditions of the 
VCO, 


n„ < 




(rad/sec) . (3-33) 

The time to reach frequency lock (not phase lock) is 


f tr "> ftcq ~ ^8FwT (w l ) (sec) (3 ' 34) 

where r = AK t\/t x and w L is the loop fiducial band- 
width (see Chapter 5). We assume m » r 2 . Typical plots 
of lock-in behavior are illustrated in Figs. II-3-3 and 
II-3-4. 
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o ■*/ z 

<t> 


Fig. 11-3-3. Leck-in behavior of a second-order loop with imperfect 
integrator, Fit) = (1 + r 2 s)/(l + ns), for Clo/AK — 0.4 and 
AKr*/n = 2. The upper bound of Eq. (3-33) is 
Clo/AK < 0.693. 
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<*> 

Fig. 11-3-4. Lock-in bohavior of a second-order loop with imperfect 
integrator, F(s) = (1 + t £ s)/(1 + ti*), for n c /AK = 0.9 and AKTf/r t 
— 2. The upper bound of Eq. (3-33) is Q n /AK < 0.693. Note that 
lock-in occurs only when the trajectory happens to pass through 
the “slot." Otherwise, the trajectory enters the periodic 
frequency lag region shown. 
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Fig. 11-3-5. Phase-plane trajectory of a second-order 
loop with perfect integrator to a doppler-rate 
Input A 0 for AK r*/ri = 2, and A 0 = AK/2r x 
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3-F. Locking the Second-Order Loop 
with Perfect Integrator 

With a loop filter of the form 


| A„ | < k(r) p 
sgn [fi(0)] = - sgn (A„) 


(3-43) 


When r = 2, the value k(r) in (3-43) is about 0.5; for 
(^•41) other values of r, factors are given in Fig. II-3-6. 


where n and r 2 are in seconds, the doppler polynomial 
d(t) = d n + 0<,f + % A 0 f 2 produces a constant steady- 
state error given by 



(3-42) 


where /? is the loop natural frequency (see 5-20). To guar- 
antee lock-in for r = 2, it is necessary to have (see Fig. 
II-3-5) 


When the VCO is swept for lock-in with 0„ ^ 0 and 
A 0 = 0, the equations above restrict the VCO input 
sweep e by 


v de, 
Kyco W \ 


< Hr) 0* 


(3-44) 


These relations are approximately true for the passive 
loop of Section 3-D when n » r 2 . 



Fig. 11-3-6. Normalized dopplor rates A<>//J 2 and A 0 /w 2 for which lock is 
guaranteed in the absence of noise, as a function of the loop 
parameter v = AKt^/t i. Note that for a fixed bandwidth w L , 
the optimum value of r is about two ({ = 0.707). 
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CHAPTER 4 

BEHAVIOR OF PHASE-LOCKED LOOPS WITH STOCHASTIC INPUTS 


With input and VCO noises n(f), n v (t), the basic equa- 
tion is 

kt) = sin [0(t) - k*)} + n(t) 

r r 

+ h£± n v (t) (4-4) 


and equivalent mathematical models are shown in Figs. 
II-4-1 and II-4-2. 

Owing to noise, the loop skips cycles at a certain rate, 
executing a nonstationary random walk between lock-in 
points. The actual phase-error variances are thus non- 
stationary, growing without bound as time goes on. When 
phase angles are reduced modulo 2tt, however, the re- 
sulting steady-state phase error process is stationary. 


INPUT 
n( t) 



t+ n v (t) 


Fig. 11-4-1 . Exact mathematical equivalent of the phase- 
locked loop. Sources of external, as well as 
internal, noises are shown; VCO tuning 
voltage e is also indicated. 


n(t) 
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Fig. 11-4-2. Equivalent exact mathematical model of phase-locked loop, 
with explicit reduction of <f> (mod 2«-l. The VCO has been replaced 
by K y oo/s (tuning bias and VCO noise omitted). 
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CHAPTER 5 

THE LINEARIZED ANALYSIS OF PHASE-LOCKED SYSTEMS 


5- A. Behavior of a Linear Loop 

When <f>(t) is very small, we can approximate sin </>=*<£. 
In such a case there is a linear equation relating the input 
phase function 0(t) to the loop estimate 6(t), 


There are two effects in 9(t): doppler d(t) and modula- 
tion <p(t). As a result, there are three kinds of phase errors 
in the loop; the mean-square phase-error 2 s is of the form 

2 2 = n*(t) + 8 2 + <r 2 (rad 2 ) (5-6) 

where n(t) is transient (or doppler) distortion due to d(t), 
8 2 is modulation distortion, and <r 2 is phase noise: 


omitting VCO noise. An equivalent circuit appears in 
Fig. II-5-1. Hence we can define a loop transfer function 


L(s) = 


AKF(s) 
s + AKF(s) 


m = 


sL(s) 

AK [1 - L(s)] 


(5-3) 


and L(s) has some fiducial bandwidth w L (see (2-37)). 


s,= -tL 


1 — L(ju>) I 2 S n (j*) dm (rad 2 ) 


(5-7) 


= (rad 2 ). 


A 2 


Total transient distortion « 2 is defined by the integral 
« 2 = j ii?(t)dt (rad-sec) 2 (5-9) 
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and the total phase-error E 2 is defined as 

E 2 = A- + 8 2 + <r 2 (rad 2 ) . (5-22) 

Here A 2 *s a Lagrange multiplier, in units of seconds* 3 , 
which can be evaluated in terms of loop bandwidth. 



Fis) 


I + r 2 5 

I +■ r,5 


WHERE r 2 = f? z C 

r, =(/?,+/? 2 ) C 


Fig. 11-5-2. Passive-integrator loop filter 


/ 


5-B. Calculation of Loop Bandwidths 

1. First-Order Loop 

IfF(s) = 1, 

w L = W,, = (cps) (5-11) 

W-jItfe- < 5 - 12 > 

2. Second-Order Loop, Passive Integrator 
With a loop filter of the form (Fig. II-5-2), 


The maximum gain of the loop occurs at <o, MX 
<w = -[(l + 2r)*-l]* 

T2 

= /,[ <! + 2 f~l ]“ (rad/sec) 

L (5-15) 

-* 0 (as r -* oo ) 

— ► ft (as r — *0) . 

The maximum value of L 2 of loop power gain is 


F(s) = 


1 + r 2 S 
1 + T| S 


L(s) = 




1 + r 2 




+ (r?/> 


(5-13) 


L 2 = | m«) | max 


r 2 

2 (1 + 2r) v * + f 2 - 2r - 2 * 

(5-16) 


These effects are shown in Fig. II-5-3. 


where r is the loop parameter ratio, 


r = 


AKrl 

Tl 


(5-14) 


We shall assume rr i » r 2 in the formulas that follow. 
The loop fiducial bandwidth is 

w L = (cps) . (5-18) 

The loop has a damping coefficient ( and natural fre- 
quency ft given by 




2r^ 

ft = — t tv L (rad/sec) . 


(5-20) 

(5-20) 


Fig. 11-5-3. Variation of maximum loop response l', 
noise bandwidth W L , fiducial bandwidth w L , and fre- 
quency at maximum loop response, for second-order 
phase-locked loop, as a function of the parameter 

r — AKr* /t, 
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3. Second-Order Loop, Perfect Integrator 

With a loop filter of the form 


F(s) = 


1 + T, .9 
T, S 


(5-21) 


the values of to,,, (, fi, w„ mx , and L 2 are the same as those 
given in Section 5-B3 when t_./t, — * 0 in such a way that 
r in (5-14) remains fixed. 


5-C. Optimization of Loop Parameters 

The Weiner optimization of the total phase error 
-• A 2 < 2 + 8 2 + it 2 is achieved by choosing 


where N„ r is the white-noise-speclrum portion of 
K\con,(t), and N t( . is the coefficient of the l/f noise com- 
ponent of K V con r (t): 

K v <-0 S«, .«„(/«>) = N,, v -f- 2irN , „/ 1 ( 1 , 1 (radVsec 2 eps) . 

(5-25) 

The function g(r) is given in Fig. II-5-4, but to good 
accuracy we can use g(r) = 1.55 for 3 < r < 10. 


Optimum values of r and tc,, that minimize <r can be 
found: a- is fairly insensitive to r, and any value in the 
range 3 < r < 10 is acceptable; can be found as 
the solution to 


^- 0 28e (^r)'“'- 3, (-i 5r)“° 


(5-29) 


L„ p ,(s) = i _ 


N'yA 

[S(*)Y 


(5-2.3) 


S(s) = A -£ [D(s)D(-s)] + S w(s) + No/A 2 . 


5-D. Effects of VCO Noise 

The total phase noise in the linear loop due to input, 
as well as to VCO noise, is 



(5-27) 



Fig. !!-5-4. Factors governing relative contributions of 
VCO noist to output phase noise 
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CHAPTER 6 

OPTIMIZED DESIGN OF TRACKING FILTERS (LINEAR ANALYSIS) 


6- A. Tracking Loop Design 

The filter L lip ,(s) that minimizes when the doppler 
polynomial d(t) has a Laplace transform D(s) given by 


D(s) = ^ + 

' S S 2 «•’ 


= ^ (6-1) 


(the decree of Q(s) is less than N) is given by the Wiener 
formula 




[(-i)v y + ( 'jfpj e i()(s)g(- S )]J ‘ 

(6-2) 


6-B. Optimum Filter for Random Phase Offset 

When d(t) = ft,, where ft, is a uniformly distributed 
random phase offset, we find 


Lnpt ($) 


XAw 


i 


XAir 

(W 7 


(6-6) 


W,. = w L = 


XAw 

2(3N„)“ 


(cps) . 


(6-7) 


6-C. Optimum Filter for Frequency and 
Random Phase Offset 

When d(t) = ft, + 0 0 f in which ft# is uniformly dis- 
tributed over (—ir, w), it follows that 



where the loop natural frequency p is given by 


P = 



( 6 - 12 ) 


This L„ t ,,(s) can be realized with F(s) = (1 + r 2 s)/r,s by 
setting 


r =L(* + l£-Y= IL_ 

* P \ 30J ) n„[3r(r - 


2 )] * 


AK 


3n„- 


— = j8- = — 2- (r ~ 2) 

T > ir- 


r - 2 + 


3n-: 


u., =. iii. n.[3,(, - 2))’* 


ft, is random. 


(6-14) 


When ft, = 0 above (not random), the value of p given 
in (6-12) remains unchanged, but the loop parameters 
are somewhat different: 


2 “ 

T ’ = T 

Tl 

r = 2 (so { = 0.707) 

36 

Wl 2 ( 2 )“ 


when ft, = 0 . 


1. Choice of Parameters 

If we design the loop optimally to lock onto an 0„ 
appearing at the bandwidth edge (i.e., 0 0 — 2 mhi), the 
proper values of r and { for ft 0 random are 


r = 2.282 
{ = 0.755. 


(6-19) 


Hence the loop configuration is described by 


t*. = 1.643/t v L 

AK iV i 

r, 1.180 


( 6 - 20 ) 


, , v 1 + 1.643/iej,» 

° pi w ” 1 + (1.643M> + (1.18/u>‘)*» ‘ 

(6-21) 
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The resulting transient error is 


v, - r 1 (r — 1) 
-r(opt) 0n„ [3r(r - 2)']’^ 


(r- 25) 


= 0.537 /tVi, for r = 2.282 and n„ = i rW/, . 


3. The Effects of Doppler Rates In Second-Order Loops 

If there is a small doppler-rate term A„(rad/sec 2 ) in d(t), 
the loop with F(s) = (1 + t..s)/t,.s has a steady-state error 


P** rate! 


A„Tt 

AK 


A„ (r + 1)» 

4 rtvl 


(rad) . 

(8-28) 


When the imperfect-integrator loop filter F(») ~ (1 + r 2 s)/ 
(1 + r,s) is used, there is a steady-state growth in phase 
error 


A A„ (r -4 1)*/ t 
^AK*- 4r uj \r, 

(6-27) 

These values are minimized, ie* - a fixed w'[ and r u by 
choosing r = 1. Parameters must be chosen to maintain 
less than about 30 deg error in the loop, and the loop 
must be periodically retuned. Representative types of 
transient behavior are illustrated in Fig. II-6-1. 

4. Comments on the Choice of r in Second-Order 
Tracking Loops 

An r between 6 and 10 minimizes the effects of VCO 
noise, whereas a value of 2.282 (or 2 if d,, — 0) minimises 
the total phase error in tracking a frequency offset, und, 
further, r — ■ 1 provides the best doppler-rate tracking 
capability. However, loop performance is generally rather 
insensitive to changing values of r. 
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(c) FIRST-ORDER LOOP, PHASE STEP 


(b) IDEAL INTEGRATOR SECONO-ORDER 
LOOP, FREQUENCY STEP 




SLOPE = 


ZERO STEADY -STATE 
ERROR 


(c) FIRST-ORDER LOOP, FREQUENCY STEP 


(d) PASSIVE INTEGRATOR SECOND-ORDER 
LOOP, FREQUENCY STEP 


SLOPE 


<P = 

AK 


(•) FIRST-ORDER LOOP, 

DOPPLER-RATE INPUT 

1 

r \-*.a 0+ v+Jv B 


C*. 

AK AK 

i i 


(d) PASSIVE INTEGRATOR SECOND-ORDER 
LOOP, DOPPLER- RATE INPUT 


t I Q . O t r\ *.l» * ^ 


e:6 0 +si 0 r 


Fig. 11*6-1. Response of first- and second-order loops to various inputs 
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CHAPTER 7 

DESIGN OF SYNCHRONOUS DETECTOR AGC SYSTEMS 


7-A. Synchronous Detector AGC Loup 

With the AGC device in Fig. II-7-1, the AGC voltage 
c(t) is given by 


7-C. Calibration Equation 

The steady-state mean value of AGC voltage c is re- 
lated to the input signal strength a by 


c(t) = 

C(p) £a(f) + K D n ;(«) + 20 log cos *(f) 


-k 8 J 

(7-8) 


a = K rec + 20 log (^j + c + (10 log ey 

(db-volts 2 ) (7-10) 


where 

a(t) — 20 log A(t) (db-volts 2 ) 

C(P) = K.i [1 + K°Ly(s)] 

v __ Ka€c,Kc 
Kagc “ 20 log c 

and Y(s) is the ACC loop filter. The other quantities are 
given in the Appendix with units. To convert a to db- 
watts, subtract 10 log R in from a, with R in in ohms. An 
equivalent circuit is shown in Fig. II-7-2. 


Thus, given c , one can infer a value of a . A typical cali- 
bration is given in Fig. II-7-3. The AGC voltage fluctu- 
ates with noise; its steady-state variance is approximately 

var [c] = £(20 log e) 2 

+ 2(10 log e)V-^Jc 2 (0). 

(7-9) 

where a 2 is the phase noise variance (5-7) . 



Fig. 11-7-1. A synchronous-detector AGC loop, using a phase-locked loop to provide a coherent reference. 
(The actual phase-locked loop may be part of the gain-controlled receiver.) 
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METER 


I 

I 


Fig. 11-7-2. Equivalent diagram of an AGC loop. Linearized analysis follows 
the assumptions that A*(f) is exponentially related to c(t) and that 
a(f)»*a*(f). The input is ait) = 20 log Ait). Adjusted 
loop gain a*(t) = 20 log IA*it) e G /K»]. 



Fig. 11-7-3. Measured AGC curve showing departure 
from linear behavior 
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CHAPTER 8 

THE DOUBLE-HETERODYNE PHASE-LOCKED RECEIVER 


8- A. Basic Configuration of the Receiver 

A block diagram is shown in Fig. II-8-1; the input and 
output phases are related by 


9 = aKdKv j sin [0 - $+ (1 + M 1 )d 1 + 0 2 ] 
+ n v (t ) . (8-6) 


This is the same form that governs the behavior of a 
simple loop when we set 

AK = aKdKycoM F (sec') 

a (8-8) 

4=0-0 + (1 + M,)8, + 0 a (rad) . 

The term F is the dc gain of F(s), viz., F(0). 

8-B. Effects of Band-Pass Limiting 

The parameter a above is the signal suppression factor 
derived by Davenport (Fig. II-8-2). In addition to this, 



* 


0 * 0 (/) = /2 rfU)eo*[to>o' + 0<'>]+"O<'> 

^ 2 ^ vW)~'TZ cos f + 0 it) j 

(3) = *ri{V2 cot [(u>o _ «>ai) + 

0 v z (t) ~ Jz cot (u hz t- fi|) 

( 3 ) /(0 s </2 A cos [(ui()- <*>*i - <*>A 2 ) t + 8 - 0 + 0 |] + it/ it) 

0 IT a cot [(u) 0 - toi A1 - W A2 ) t + 6-9 + 0,J 
+ «<<'>} 


0 y it) s Kd {a tin [(w 0 - - (I + * 1 ) <»/,z) t + 6-8 

+ (l+ M\) A| + ®z] + » </>} 

0 * (/) s Fip) y it) 

0MM~V2co.[^+ /(/)] 

0 * 3 (/)~ V 2 tin(V|u A2 t + Af,9| - 6 Z ) 

0 *A6C<'> s *«tf + co * [("0~ "*l ~P + + ®* 

+ (I + 4/, ) 6, + 0 2 J + n m} 

C ( / ) s K( Yip) /yyGC i 1 1 


Fig. 11-8-1. The double-heterodyne receiver, with equations for signals at each point 
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Fig. 11-8-2. Davenport’s band-pass limiter zonal SNR 

curve 

the ratio A-/N 0 at the limiter input comes out of the limiter 
as A-/N„V, where F is a factor shown in Fig. II-8-3. 
With the loop filter F(s) = F ♦ (1 + + r,s), we 

thus have (linear theoretic approximation): 


Fig. 11-8-3. The ratio of band-pass limiter output signal 
to-noise spectral density to that at input. The 
reciprocal of this curve gives the limiter 
performance factor T 


1 + 0.345p* 
0.862 + 0.690p» 


(8-18) 


, _ n„w l r 

a — — 7- — 1 = 


■ - m(l + f .)'( 1 + -^ r ') < rad! > 

(8-19) 


aKdKvcoAf Ft; a 

r = - — r 0 

ti a 0 


(8-14) 


to,, = 


1 + r _ 

2r.» Wl ° 


1+ (i) r - 


1 + r 0 J 


(cps) 


(8-15) 


r* _ /_a_V 

2 v «o ; 


(5-20) 


kT 


No = -g- (watts/cps) 


(2-65) 


Al = N 0 w Ln (watts) 



_ ( 0.7854 p H + 0.4768p * 

“ _ \ 1 + 1.024p„ + 0.4768p j, ) 

,, _ ^(KdKycoM F)’tI _ r 2 0 

8 r\w H r 0 + 1 


( 10 - 1 ) 

( 10 - 2 ) 

(10-3) 

(8-13) 

(10-5) 

(10-4) 
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CHAPTER 9 

BEHAVIOR OF PHASE-LOCKED LOOPS (NONLINEAR ANALYSIS) 


9- A. The Spectral Equation 

The phase-noise spectrum can be approximated by 


S**(s) — 


K 2 N 0 F(s)F(-s) 

—s 2 +i)AK [sF(-s) - sF($)] + (AKy) 2 F(s)F(-s) 

(rad 2 /cps) (9-11) 


v 


e -o 2 /2 


(9-14) 


y = 



(9-19) 


when o- 2 < 1 (rad 2 ). Then (9-11) can be integrated to 
give a 2 : 


N O^L(eq) 

A 2 y 2 


(rad 2 ) 


(9-28) 


where Wu^) is the loop equivalent bandwidth, 




2b- 



S** (/<!)) 


(cps) . 
(9-27) 


9-D. Behavior of First-Order Loop 

For F(s) — 1, a better value of y than (9-19) is (see 
Fig. II-9-2) 



for a 2 < 1. The phase-noise variance is thus a transcend- 
ental equation in <r-‘, whose results are given in Fig. II-9-3, 

<r e-° 2 ' 2 (sinh <r)* = ■ (9-25) 

At tr = 1, the value of N„ vcJA 2 is 0.857. 




0 1.0 2.0 3.0 4.0 


N q w l /A Z . NqK/ZA 

Fig. 11-9-2. Variation of tho paramotor y as a function of 

the Gaussian variance a 2 , for various forms of ft** (r). Fig. 11-9-3. Comparison of linear, quasi-linear, and 

Note that, for a 2 < 1, there is not a significant linear-spectral approximate methods with the 

dependence on the form of ft**(r). actual behavior of the first-order loop 


1 
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LINEAR THEORY 


NONLINEAR THEORY 


■ 

I 

I 


■ 

I 

I 


SIGNAL STRENGTH, db ABOVE A% * N 0 w Lo 

Fig. 11-9-4. Comparison of linoar and nonlinoar thoorios for socond-ordor, constant linoar bandwidth loop, i. 

♦ho valuo of r is kopt constant at r 0 = 2. 






PHASE NOISE VARIANCE <r z rad 
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Fig. 11-9-5. Comparison of phaso-noiso variancos by iinoar and linoar-spoctral approximations. Tho noiso density 
it fixed, and the signal level it varied. The value of r it taken as r 0 — 2 at a reference tignal level of 
A* = 3N 0 /2t 2 = N a Wi, n . Note that the ultimate roll-off it 5db/decade, rather than lOdb/decade, 
as in Fig. 9-4. Note alto that even the linear approximation producet some curvature of near AJ. 
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9-E. Calculation of Behavior of the 
Second-Order Loop 

With the usual passive-integrator filter second-order 
loop, 


WLjeq} _ 1 + ry 

(1 + ,) [i + i!^]" 

t-4 l+ VT- 


(9-29) 


(9-33) 


These, with (9-28) above, specify the loop behavior. 



Figures II-9-4, II-9-5, and II-9-6 indicate the true loop 
behavior compared with that predicted by linear theory. 


Fig. 11-9-6. Variation in bandwidth and damping param 
eters as a function of signal strength. The value of r 
at a reference signal level AJ* = 3N„/2 t 2 = 

N 0 w L() was taken as r 0 = 2. 
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CHAPTER 10 

DESIGNING A DOUBLE-HETERODYNE TRACKING LOOP 


70-4. Definition of Receiver Threshold 

Threshold is defined as a signal level A 2 0 such that 18 

A* = N n wr, a = — 2 T -f-- (volts*) (10-1) 

(assuming r 0 Ti » r 2 ). The loop is not acting linearly at 
this signal level; nevertheless t v,.„ is the bandwidth a 
linear loop would have at A = A„. 


70-8. Tracking Loop Performance of the 
Double-Heterodyne Receiver 

The ratio (A/A 0 )- is the receiver margin: 

m = (A/A 0 ) 2 . (10-2) 

With measured parameters Wn,KjK y C oMF, r 2 , and n, the 
value of r 0 can be calculated (when 10tt>£, 0 < u>«) as 


ir(KdKvcoMF) 2 rl 
St- W„ 


(10-4) 


r o + 1 


y = (1 - e- a? )/<r 2 
v = e- a 'U- 


w r , = 


1 + r 

~2 WT = «*. 


i _i_ n 

1 H r„ 

(Xo 

1 + r n 


(10-6f) 

(10-6g) 

(cps) (10-6h) 


t^L(eq) ~ 


1 + Ty 


= 


4 i+j ^r 


1 + r o 


(10-6i) 




N oti^L(cq) ^ ^ ^ 

A 2 y 2 my 2 


1 + 


1 + 1*0 


(rad 2 ) 

(10-6m) 


B = = (-i-Vi^L (10-fln) 

P r + 1 u/ r„ + l 


The only modifications to the linear theory in Chapter 8 
that need concern us for a- < 1 are that A is replaced by 
Ay, r by ry, and w L by tc L(eq) . The resulting equations 
that specify loop behavior are 


p H = m 



a 


0.7854p„ + 0.4768p 2 “| * 

1 + 1.024p„ + 0.4768p 2 


1 + 0.345p* 
0.862 + 0.690p„ 


(10-6a) 

(10-6b) 

(10-6c) 


aKjKycoAfFr 2 

T» 



(10-6d) 


“Since N» = kTR/2 (voltsVcps), the threshold signal power is 
P« = AJ/R = kTw,./2 = kTbi. (watts). 



r 0 , obs 

Fig. 11-10-1. Variation in margin producing a 2 = 1 as a 
function of threshold design parameter, 

to ® BuitjKvcoWfTj/Ti 
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7 0-D. The Signal Level Producing <r 2 = 7 

The signal level A~, producing a 2 ~ 1, corresponds to 
a margin m, approximately equal to 


m, 


( 


r 0 r ] 

2y,(ro + 1) 



1 + 


4 (f, + 1) 



2 


3.13 (for r„ = 2). 


( 10 - 8 ) 


Hence the <r = 1 margin is about 5 db at r„ = 2. (This 
result is shown in Fig. II-10-1.) 


10-E. Choice of Receiver Parameters 

Choosing values of N 0 , w Lo ,K,i, K V co, MF, and u>„ 
specifies a typical receiver design. However, the only 
parameters required to plot performance are m, y,„ pn 0 
and perhaps t 2 /t, \tJt x is r.ol needed if r, » r-..). Such 
plots are given in Figs II-10-2 and II-10-3, 



Fig. 11-10-2. Comparison of linoar and nonlinear approximations to loop rms phaso orror, as a function 

of loop margin 
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MARGIN m, Ob 

Fig. 11-10-3. Variation of loop bandwidths and damping factors, os a function of loop margin 
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APPENDIX 

Nomenclature 

Those symbols used throughout the text are listed here, along with assigned 
names, proper units, and the equation number that either defines the quantity, 
gives its first use, or else refers to its approximate point of first introduction. 
Quantities labeled “abs” are dimensionless. 


a either mean value of a(t) below, or square 
root of a 2 below. 

a* coefficients in expansion of p(<j> u (9-5). 

a 2 variance of <!>(<), in rad 2 , (9-4). 

a(t) input rms signal A(t) expressed in db-volts 2 , 
(7-4). 

di(t) that part of a(t) due to deterministic and non- 
stationary changes in signal, in db-volts 2 , 
(7-11). 

av(t) iha‘ part of a(t) due to stationary random 
signal fluctuations, in db, (7-11). 

cZt(s) the Laplace transform of a d (t), (7-13). 

a*(t) adjusted receiver attenuation A*(t) expressed 
in db-volts 2 , (7-4). 

A rms signal amplitude into loop, volts, (3-1). 

A(t) rms signal amplitude into receiver at time f, 
in volts, (7-1). 

A*(t) receiver attenuation factor at time t, abs, (7-1). 

A,i(t) nonstationary part of A(f), in rms volts, (7-12). 

A,/,(t) stationary gain fluctuation part of A(t), abs, 
(7-12). 

A t receiver rms signal level producing a = 1 
rad, volts, (10-7). 

b H fiducial (one-sided) bandwidth of the linear 
transfer function H(s), cps, (2-37). 

B h equivalent noise (one-sided) bandwidth of 
the linear transfer function H(s), cps, (2-27). 

c, c(t) AGC control signal, volts, (7-3). 

C capacitance, farads. 

C(s) AGC closed-loop transfer function, (7-8). 

90 


/ 

d(t) doppler phase function on input signal, rad, 
(3-11). 

di(t) time-varying part of doppler phase function 
on input signal, rad, (5-9). 

D(s) doppler phase function in s-domain, D(s) = 
_£»)], (3-12). 

e VCO tuning bias, volts, (3-5). 

e a AGC gain-adjust bias, volts, (7-1). 

E( ) statistical expectation operator (2-1). 

f frequency variable, cps. 1 9 

ft first IF frequency, cps, (8-1). 19 

/ 2 second IF frequency, cps, (8-1). 19 

f hl first-mixer heterodyne frequency, M times 
VCO output frequency, cps, (8-1). 19 

fh 2 second-mixer heterodyne frequency, from in- 
ternal oscillator, cps, (8-1). 19 

F (finite) dc gain of loop filter, abs, (8-8). 

F(s) linear loop-filter transfer function, (3-7). 

3 Fourier transformation operator. 

g conductance of diode, or transconductance of 
triode, mhos, (2-66). 

g(r) VCO noise term coefficients, abs, (5-28), 

G conductance, mhos, Fig. 2-9. 

G xr (/o)) single-sided spectral density of x(t), volts 2 / 
cps, (2-26). 

h Planck’s constant, 6.625 X 10 _S4 joule-sec, 
(2-53). 

h(t) un’t-impulse response of linear filter H, volts, 
(2-14). 

’’May take on a negative value. 
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H a filter operator. Section 2-B. 

H(s) linear filter transfer function, X[h(f)l> (2-26). 

i(t) current, amperes, (2-54). 

/ V — 1 , the imaginary unit. 

' k Boltzmann’s constant, 1.38 X 10 -23 joule/ °K, 

(2-53). 

k(r) maximum normalized lock-on doppler rate, 
abs, (3-43). 

K K,K m K V co, equivalent simple-loop gain, volts ’sec -1 , 
(3-8). 

Ki rms VCO signal output, volts, (3-2). 

K a receiver AGC attenuation, db/volt, (7-3). 

K AC c AGC-loop equivalent loop gain, abs, (7-6). 

Kc AGC amplifier gain, abs, (7-6). 

K* phase detector gain, volts/rad, (8-2). 

K 0 AGC detector gain, (volts peak out)/(volts 
rms in), (7-3). 

K n first IF gain, abs, (8-1). 

K /2 second IF gain, abs, (8-1). 

Km mixer gain (simple loop), volts' 1 , (3-4). 

Kg adjusted receiver attenuation with no AGC, 
db, (7-3). 

Ky co VCO gain constant, rad/sec-volt, (3-5). 

K r . e actual receiver attenuation with no AGC, db, 

(7-3). 

I limiter rms output level, volts, (8-2). 

L inductance, henries, Fig. 2-9. 

L 2 maximum value of | L(/w) | 2 , abs, (5-4). 

L(s) linear loop phase-transfer function, (5-3). 

X Laplace transformation operator. 

m receiver margin above threshold, abs, (10-2). 

M VCO output frequency multiplication factor, 
abs, (8-1). 

Mi internal oscillator frequency multiplication 
ratio, (fi/fh >) ~ 1, abs, (8-4). It must be a ra- 
tional number, and may be negative. 

n 0 (t) receiver input noise waveform, volts, Sec- 
tion 8-A. 

n 4 (t) loop input noise waveform, volts, (4-1). 


n(t) simple loop baseband noise waveform, volts, 
(4-3). 

N r (t) equivalent VCO noise referred to VCO input, 
volts, (5-24). 

N noise power out of linear filter, volts 2 , (2-23). 

N 0 noise (two-sided) spectral density, volts 2 /cps, 
(2-20). 

N 0 o receiver input noise (two-sided) spectral den- 
sity of n 0 (t), volts 2 /cps, Section 8-A. 

N. noise (one-sided) spectral density, volts 2 /cps, 
(2-27). 

N or VCO adjusted input white-noise density, 
(rad/sec) 2 /cps, (5-25). 

N lv VCO adjusted input l/f noise density, 
(rad/sec) 3 /cps, (5-25). 

p Heaviside operator, d/dt, (3-6). 
p(A) probability density on A, abs, (2-1). 
p(s) denominator of F(s), (3-15). 

P input signal power, A 2 , volts 2 , (2-29). 

P y power in y(t), volts 2 , (2-19). 
q(s) numerator of F(s), (3-15). 

Qo resonance quality factor, <■ > 0 L/R, abs, (2-57). 
Q(s) numerator polynomial of D(s), (6-1). 

r second-order loop-parameter ratio AKt 2 / t,, 
abs, (5-14). 

Rx resistance, ohms, Fig. 2-9. 

Rx V (ti, t-i) statistical cross-correlation between x(t t ) and 
y(t 2 ), volts 2 , (2-4). 

R„(r) statistical autocorrelation of stationary proc- 
ess x(t), volts 2 , (2-5). 

3lu(r) time autocorrelation of the function x(t), 
volts 2 , (2-6). 

s complex frequency variable, rad/sec. 

S output signal power of linear filter, volts 2 , 
(2-33). 

S„(/<») two-sided spectral density of stationary proc- 
ess x(t); £F[iWr)], volts 2 /cps, (2-8). 

S„(b) two-sided spectral density of function x(t); 
£?[iR**(T)], volts 2 /cps, (2-9). 
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S(s) spectral density of total phase-error, rad 2 /cps, 
(5-23). 

S'(s) spectral density of total AGC error, (db- 
voltS')'-/cps, (7-13). 

t time, sec. 

t u f 2 specific instants of time, sec. 

f acQ first-order loop phase-acquisition time, sec, 
(3-27) 

f f r eq acq second-order loop frequency-acquisition time, 
sec, (3-34). 

T,T n equivalent noise temperature of receiver, °K, 
(2-53). 

T c cathode temperature of diode or triode, °K, 
( 2 - 66 ). 

U receiver measured-parameter ratio, abs, ( 10-4) . 

v(t) VCO output, volts, (3-2). 

Wh two-sided fiducial bandwidth of linear filter 
H(s), cps, (2-37). 

w L two-sided fiducial linear loop bandwidth of 
phase transfer function L(s), cps, (5-3). 

u>k«d equivalent two-sided fiducial bandwidth of 
loop output spectrum, cps, (9-29). 

Wh two-sided noise bandwidth of a linear filter 
H(s), cps, (2-24). 

W L two-sided noise bandwidth of a linear loop 
whose phase transfer function is L(s), cps, 

(5-4). 

x(t) arbitrary time function. 

X(s ) Laplace transform of x(t), X [x(f ) ] . 

y(t) time function. 

Y (s) AGC loop-filter response, (7-5). 

z(t) time function. 

Z(/<o) complex impedance function, ohms. 

a signal voltage suppression factor, abs, (8-1). 

a« value of a at threshold, (10-6). 

/? second-order loop natural frequency 
(AK/t,)'*, rad/sec, (6-12). 

y(s) sinusoidal-error spectral ratio function, abs, 
(9-7). 

v y(0). 


y tp sinusoidal-error function total-power transfer 
ratio, abs, (9-15). 

8 2 modulation distortion, rad 2 , (5-6). 

Siock first-order loop “in-lock” constant, rad, (3-27). 

S T (t) rectangular unit-pulse function, volts, (2-13). 

8(f) Dirac delta function, lim 8 T (f)asf-» oo,(2-13). 

tij, total transient distortion, rad 2 , (5-9). 
loop tracking error, rad, (8-7). 

£ linear loop damping factor, abs, (5-20). 

£ (e q) nonlinear loop equivalent damping factor, 
abs, (9-33). 

tj Bussgang coefficient, abs, (9-3). 

&i phase angle, rad. 

$„ 0(0), initial value of loop phase offset, rad, 

3-13). 

0(t) or 0 input signal phase function, rad, (3-1). 
t ) or 6 loop estimate of 0(f), rad, (3-2). 

A 2 Lagrange multiplier, sec -1 , (5-22). 

An d ,n> (0 +), the nth doppler moment, rad/ 
sec", (3-15). 

fi mean of a stationary random variable, volts, 

( 2 - 12 ). 

n(t) mean of a nonstationary random variable, 

^ volts, (2-10). 

*• 3.14169... 

p v signal-to-noise ratio of y(t), abs, (2-35). 

<T= variance of stationary random variable x, (x 
units) 2 , (2-12). 

as (f) variance of nonstationary random process x(t) 
at time f, (x units) 2 , (2-11). 

<r 2 variance of loop phase-noise, rad 2 , (5-7). 

<72 variance of loop phase-noise due to noise in 

VCO 

VCO, rad 2 , (5-26). 

r f! — f 2 , time difference, variable of R«*(r), sec, 
(2-4). 

ri second-order loop denominator time constant 
of F(s), sec, (3-28). 

r 2 second-order loop numerator time constant of 
F(s), sec, (3-28). 
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<f>(t) = <f> two-heterodyne loop detector error (same as 
for simple loop), rad, (5-8), (8-6). 

<j>$$ steady-state loop detector phase error, rad, 
(3-12). 

</> o initial value of loop detection phase error, rad, 
(3-32) 

$ stationary equivalent phase-error process, 
rad. Section 4-B. 

<0 angular frequency variable, rad/sec. 
o) 0 VCO short-circuit output frequency, rad/sec, 
(3-1). 

“max frequency at which |L(/u>)| is maximum for 
second-order loops, rad/sec, (5-15). 
o>x first IF frequency, rad/sec, (8-l). ,!t 
o >2 second IF frequency, rad/sec, (8-1). 19 
whi first mixer heterodyne frequency, M times 
VCO output frequency, rad/sec, (8-1 ). 1:1 


<i>h 2 second mixer heterodyne frequency, from 
internal oscillator, rad/sec, (8-1). 1! ’ 

T limiter performance factor, wip t /w H pH, abs, 
(8-17). 

A 0 initial value of doppler phase-rate, rad/sec 2 , 
(3-13). 

S 2 mean-square phase error, rad 2 , (5-6). 

2 2 total mean-square phase error, rad 2 , (5-22). 
i p(t) modulation phase function, rad, (3-11). 

O = Cl(t) 4>(t), loop frequency error, rad/sec, (3-20). 

fi 0 initial value of loop frequency offset rad/sec, 
(3-12). 

Q m ax maximum value of Cl for which loop locks in 
absence of noise, rad/sec, (3-23). 

'"May take on a negative value. 
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